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I. INTRODUCTION 


There are a number of methods in use today that numeri- 
cally solve for the currents on simple wire antenna struc- 
tures. The majority of these methods utilize an integral 
equation solution to Maxwell's Equations which also match the 
boundary conditions on the structure and satisfy the radiation 
conditions for the radiated fields. One of the most widely 
used implementations of the integral equation method is the 
Numerical Electromagnetics Code (NEC) {Ref 6]. This program 
has been proven to be highly accurate for antenna structures 
that can be modeled using individual wire segments or surface 
patch elements. 

A new method is presented here which does not rely upon 
an integral equation formulation. This technique encloses the 
antenna structure between two parallel ground planes, as an 
approximation to the free space condition, to allow a partial 
field representation using a periodic Fourier series. In the 
initial investigation being reported here, the simple monopole 
radiator is being considered. In such a case, cylindrical 
regions are established around the structure and the harmonic 
expansions of the magnetic and electric fields are matched 


along the regional boundaries to solve for the unknown 


expansion coefficients. These coefficients are then used to 
compute the current imposed on the surface of the antenna. 
This project began as a search for the viability of 
computing the currents on a top loaded monopole between two 
ground planes using cylindrical harmonic expansions in three 
regions. After a number of weeks of inconclusive results, it 
was decided to look at the simpler case of the monopole 
antenna between two ground planes as presented in Chapter 2. 
During the validation of the code, as presented in 
Chapter 3, it was found that the number of terms used in each 
region can be reduced by as much as 90% and still yield 
consistent results for the surface current and input im- 
pedance. However, the representation of the electric field 
along the surface of the antenna is highly inaccurate when the 
number of terms are insufficient to adequately represent the 
field. This phenomenon will be discussed in Chapter 4. 
Additional research revealed that a simpler method could be 
used that would remove the requirement for one complete 
expansion. This is presented in Chapter 5. Finally, con- 
clusions and recommendations are presented in the last 
chapter. In addition, Appendices are used to contain much of 
the detailed descriptions of the computer algorithms and 


resultant programming for this effort. 


A. METHOD USED 


The general approach is to place the antenna structure 
being analyzed between two perfectly conducting ground planes 
and divide the surrounding area into a number of cylindrical 
regions, using a cylindrical coordinate system, that coincide 
with the natural boundary interfaces of the antenna as shown 
in Figure 1. The upper ground plane is introduced to com- 
pletely close the regions and produces multiple images, each 
identical to the previous image but with opposite polarity. 

Once the regions are defined, we start by obtaining the 


solution to the complex scalar Helmholtz equation. Then, by 


Ground Plane 


Ground Plane 





Figure 1 The Unloaded Monopole Antenna Structure 


using an harmonic expansion of the resulting scalar wave 
function, ¥, we write an expression for the vector electric 
and magnetic fields in each region. The tangential component 
of the electric field, as approximated by the truncated 
expansion in Region I, is set equal to the "known" field in 
the gap and set equal to zero along the cylindrical surface 
of the perfectly conducting monopole. The tangential electric 
and magnetic field expansions are then matched across the 
boundary interface between Regions I and II. The resulting 
set of equations is reduced to a system of one equation with 
one unknown by sifting out some of the coefficients using the 
orthogonality principle of Fourier moment integrations. The 
resulting system is then solved for the remaining set of 
unknown coefficients. These coefficients are then used to 
compute the current distribution along the surface of the 
antenna. To validate the source code, the results are then 
compared to those obtained from the Numerical Electromagnetics 


Code (NEC) to ascertain computational accuracy. 


B. ASSUMPTIONS MADE 


In order to simplify the derivations, assumptions were 
made regarding the electric characteristics along the ground 
planes and the surface of the antenna structure, and the 
electric field produced by a driving source of constant 


frequency. First, the ground planes are modeled as perfectly 


conducting ground planes of infinite dimension. In reality, 
a ground plane can be considered to have infinite dimension, 
for purposes of antenna current calculations, if it is larger 
in dimension than about 10 times the size of the antenna [Ref. 
1}. Additionally, the electric field produced by the driving 
source is considered to be of constant magnitude within the 
gap region. This is a standard source model employed in 


integral equation modeling of antenna structures, [Ref. 1]. 


ll. COMPUTATIONAL PROCEDURE 


This chapter will present the mathematical background 
required to solve for the currents and the input impedance on 
a simple monopole antenna structure using cylindrical harmonic 
expansions. The final solution will be reduced to a system 
of one set of unknowns which will be solved using a digital 
computer system. An overview of the software code used, 
including error detection, hardware requirements and a brief 
discussion of the input and output file structures will also 
be presented. Details of the algorithm implementation are 


relegated to Appendices. 


A. THEORETICAL TECHNIQUE 
The antenna structure shown in Figure 1 (page 3) is 
orientated within the cylindrical coordinate system as shown 
in Figure 2. We begin this development with the complex 
scalar wave equation, known as the Helmholtz equation, given 
in cylindrical coordinates as [Refs. 2, 3] 
2 
V v (pz) +k y(p.z)=0 
(Eq. 1-1) 
where ~ is used to denote the electric type Hertz potential 
and not the magnetic flux [Ref. 3]. Expanding this equation, 
and noting that the structure and the fields are assumed to 


be axisymmetric, one obtains the following: 


Ground Plane 





Figure 2 Cylindrical Coordinate System 


{Eq. 1-2} 
Utilizing axisymmetry, the transverse magnetic (TM) 


solutions are given in cylindrical coordinates as [Ref. 3] 





Z 
aye. H, =0 
y Op dz 
0 
BG Hg=- 
dp 
E, =x|+kg VY H, =0 
y \a2 
{Bees 1-3} 
where 
y =/ OE, 
PEG, di 4) 


The following boundary conditionons must be satisfied: 


*- The electric field tangential to the ground planes must 
be zero (E, = 0). 


* The electric field tangential to the top of the antenna 
structure must be Zero. 


¢ The tangential electric field along the antenna structure 
must be zero (E, = 0). 


*- The electric field along the gap distance d must be a 
constant value (E, = -V,/d). 


+ The electric and magnetic fields must be continuous 
across the boundary between the two regions above the 
antenna structures (E,, Hj). 

* The radiated fields in Region I must be purely outbound 
at large distances from the antenna and must satisfy the 
radiation condition that E, = n)H,, where mn) = 3779. 
Using separation of variables, where the product solution 


has the form y~ = R(p)Z(z), one obtains the following equa- 


tions: [Ref. 3] 


dp 


ea [p £8.) [(ey0)-n" =0 


Oa 7 =() 
i (Eq. 1-5) 
The first equation is recognized as Bessel's equation of order 


n. The solutions to Bessel's equation are represented by any 


of several special functions, including [Ref. 3] 


SECS Cy eae. (Uy pemren (Ra) (Eq. 1-6) 


where J,(k,p) and Y,(k,p) represent Bessel functions of the 
first and second kind of order n, and H,"?(k,p) and H,‘”)(k,») 
are Hankel functions of the first and second kind of order n, 
where (2) = J. + 7 Yee As any two of these functions are 
linearly independent of each other, the solution to Bessel's 
equation can be represented as a linear combination of any two 
of these functions. The solution to the Z equation is a 
linear combinations of harmonic functions of the form e?*? and 


a , 


Using these results, the desired solutions to the Helmho- 


ltz equation becomes 


Jo(kp Pp) 
i, ee 


He ep 0) ' ) cos (k, z) 
Choe 1-7) 
where 

2 2 2 

ki +k, =k, 

ie 

= (ape 

(EqG@e1=8) 


and where k, is defined as the wave number, qm the angular 
frequency, # the permeability of a vacuum (4m x 10” H/m) and 
« the permittivity of a vacuum (8.854 x 10°” F/m). Acceptable 
quantized values of k, are found by substituting the product 
solution into Equation 1-3 and setting E, = 0 on the upper 


ground planes in Region I while enforcing E, = 0 on the upper 


ground plane and the top of the antenna in Region II. This 


gives, respectively, 


k, = Vian Reowvonel {Eq. 1-9.a} 
n 
l 
meee) oF = 
k, = in Region II (Eq. d-S 2b) 
2 l-h 


It should be noted here that the trigonometric terms used in 
Equation 1-7 will uniquely satisfy the first three boundary 
conditions listed previously for each term in the series. 
For large arguments, the Bessel and Hankel functions 
listed above closely resemble harmonic functions. As seen in 
Table I (Ref. 3], J,(k,e) and Y,(k,o) are analogous to cosine 
and sine functions respectively, while H,‘”(k,p) and H,‘” (k,p) 
are related to complex exponential functions. Therefore, 
Bessel functions can be used to represent standing waves, 
while the Hankel functions describe travelling waves; specifi- 


cally H,°’ represents inbound waves while H,/‘”’ 


represents 
outbound waves. Using these properties, the potential field 
in Region I can be written as the product of Hankel and cosine 
functions where the cosine term implicitly satisfies the 
tangential electric and magnetic field boundary conditions. 

Because the higher spectral bandwidth of the potential 
field generated by the gap voltage was expected to require a 


larger number of terms in the Fourier series expansion than 
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PROPERTIES OF SOLUTIONS TO BESSEL'S EQUATION 


Table I 


°O xf U ‘uO PUB wo PUY ‘9 = u ‘ed 30] pay T are suORDUD] [estag OF3 0 — y TOT AA + 


fm (OD)*T =m (S¥f)*T SMO ITN} oT} “Of — — XY TOM « 


Aovelesstopi— erent) © = 2) ye = C1Nre Pree TS 








Sone SUIPUBIS Paz{[RoOf—To] Amd F o<u o(97)e 
rc . a The 1(T — ¥)eZ . 
yer 9} 1 /$ OS oye Pees. Coney 5 pea Be 
Sppyg 3uaosemEse Oaj—AIBUISUUIT ¥ Buofe ( t yee — x" z] r (97)*N 
c Oo@= oY zoOtmnu 
eyrapuy 


ese BZUTPURIS—][Cel 
a 
aaua ZUPUTIT pazyeooj—xajdmoo 7 
[Oxia E + (7X) HI7 (°x)"L 


splay aceowetrEss OMI—ATBUTBUUT Y | or oy 


asBa SUTPUTIS—JEel 
Deen EERIE mama ee, Se hp 
(99) ee 














SuyesBI} parsnuaz20—xe]cmI0D ¥ o<u 
of — oY Y (1 — ues # (9%) 
WE 
of —— earl — (ox)™ o w 
pres gusoseuese—AISULIVUI! ¥ alm — Ox ae bz A (62) "N! (oner | 4a" A 
Od mr 
o- oy Oo u = Se) ot =F + T 
SAR. BOTSARI-PIVMINO—][el F i 
a To  ———————— 
esate 
Zuyssaey persnuez3e—redaoo 7 (CDE) eglese 
al om Sy me aye OS TS wee _ sen oy)" NE + (On)*L | (CHT 
Pres guaosauwse—AIVUUEUT! ¥ ool — OX at LZ= = (C9) N: UL 47) 
o= oy oe ay — T 














guonsiusssidaz 
eanstrez[V 


(= <— joy}) eepnti0} (ox)*¢ 
quauIngse-a2Ie J qusmIngsre-Teus 





SARK BUTIABI-PIBMUI-—[Ssl ¥ 
Leen eee Ean ees asain ——— ne en eenmeanmaal 
— (0 — 6H) vernm0} | 





uonwzasdreimy Pwormsy I gen ruyuy e027 





11 


that of the quasi-singular field at the end of the antenna, 
the total field in Region I was separated into two parts as 


shown below 














v(p.z)= >. a, H Soyo yoo om Hy W,,p ) cos aoe 
n=0 
=D 
{Eq.1-10} 
where 
2 
Da IE 
Yn = e3-(*F) 
(Eq. 1-11) 


In Region II, the potential field is constrained as a 


standing wave and can be expressed as 


y(p,z)= > C, J(u, Pp) cos oe a 


a {Eq. 1-12} 


where 


(Eq. 1-13) 
and gq=I/1-h. 


Using Equation 1-3, the electric field component at the end 


of the antenna is approximated by the truncated expansion 


12 


Ny 


Z 
” y Var HP ¢vga) cos|" “ | 
a= 6 





a (aag) — 
aM ) 


0, for OSzsh 
] 2) n 1\(z - h) y 
=) U, Cn Jg(Upn@) cos a forhs<zsl 
* n=0 
{Eq. 1-14) 


and the electric field due to the gap voltage is given by 





My 
1 2 (2 mz 
E,,(4z) ==> > Vv bm He (wv, a) cos ( 
m =0 
V o 

ie ee < 

; re for 0O<z <d 

0, for d<z<l 


(Eq. 1-15) 


where N, represents the number of terms to be used to in 


Region I, N, for Region II, and M, for the electric field 


component due to the gap voltage. 

The last boundary condition to enforce is that the 
magnetic field be continuous across the interface of Region 
I and Region II. Again using Equations 1-8, 1-10 and 1-12, 
and setting H, (az) = H,"(az), an expression for the phi 


component of the total magnetic field can be written as 
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nNMZ 
H g (4,2) = Dae sa) 0 | j | 





My 
mnz 
+ > Vn b nH vga) 008 | | 
m =0 
Ny 


-h 
me 2 UnCyd 1 (U,a) COS =] 


ny {Eq. 1-16) 
Equation 1-16 can be reduced to one set of unknowns by solving 
for the coefficients b, analytically and using this result to 
solve for the c, in terms of 6, and the unknown a,'s. This is 
accomplished by sifting out the desired coefficient using the 
orthogonal property of Fourier moment integrations which can 


be defined over the regions where the fields exist as 











| 0, me#n 
(1) _ IT 2 nT 2 | 2 
I, = Cos cos |— ge 7 oa 
0 l, m=n =0 


(Eq. 1-17.a) 


| : m #=Nn 


0 
o) : s 
ae gos(™% CA) )ogg(2# CEH) y a > Mm =H ea 


h q, m=n =0 


Pong 


(Eq. 1-17.b) 
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d d, 


a eo, 
hl = cos ( 7 j OU) = 


foonian? 10d 
y sia ) m # 0 
mim 


{Eq. 1-17.c) 








Applying the integrations to both sides of Equation 1-15, the 
b, coefficients can be solved for in terms of known functions 


resulting in an expression for the terms in the second 


summation of Equation 1-16, as shown below: 
M, 


2 
pe Ve on cos| 


nNz 
S 





dz 











m=0 











ef cos (2 
7 S 
(Eq. 1-18) 
Therefore, 
2 ey 
v, b, Hew, aI = > 1,” 
(Eq. 1-19) 
which can be rearranged to give 
(3) (2) 
WE, V1 HH, (W,a 
vb, HOW, Ss aa ae ae 


By applying the Fourier moment integrations to Equation 
1-16 then substituting in the expression for b, (Equation 


1-20), the c, term can be written as 
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Ny 


J g(u,a) 2 
“— = V,a, Hy WO) Ten 


2 
Uy, Gao CU) BE 
le J U,4@) k=0 


M 3 2 
Po Vie nic era a) 
- je, ) —~—- —-—_T;, , 


] 2 
ome If Ho \v,a) 








(Eq. 1-21) 
where 
I 
Cy ae cos( Fos] Mc -h | 
q 
h 
] 
a | l ees 
= ; 5; sin 
= ea 
a 
(Eq. 1-22) 


We now have one equation with two unknowns. A second expres- 
Sion for the unknown a,'s can be obtained by applying the 


moment integrations to the Equation 1-14 such that 


N, 
2 
» Vv a,Hy (V,a) | cos we Cos ee dz 


n=0 0 


N, pl 
= » wu CyrJo(u,a) cos = 2) cos ("dr 


n=0Q 


h 
(Eq. 1-23) 
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which gives 


N 
2} 2 
(2) 1 2 
Vn Ging Vint Im = SY) Un Cao Und) Tiny 
n=Q0 
(Eq. 1-24) 
Rewriting Equation 1-21 yields 
N, 
2 2 
UnCnI (Un Vn = >) Vy ae (Vy a) Tin 
k=0 
Vg et te) Hv, a) 
5 i). (yc 
Sv tl? H,?%,¢) 
(Eq. 1-25) 


Substituting this result into Equation 1-24 yields 


Ny 


eat ys Q) 
n.d 
VedmH OW, ate = Yu, V4 Hy Va) Ten 
ee; fs aaa) k =0 


M | 3 (2) 
ae Che) nn Jae Poe EL One 
“JO Eq “a me ae 
17 Wa) 1c a) k= =0 Yk th: 1 H 9 va ) 
{Eq...1=26} 


Collecting terms and rewriting yields 


Ite 


N 


N 2 
eRe eel (9 1 if ey be 
> vy a H\(v, a ) ee Vm ia Aen. ay ie 
k=0 n=0 Jn Jy Uy) 


M 1,2 (2) N 
= j@ io pi al (V4 ) ; Undo (u,,@) | aap ye: 
aa0ae 1, PH 


Ce) | an Re Ge 


(Eq. 1-27} 
This can be written in matrix form as 
N 
» Am4ae = Br for m =0,N, 
eal (Eq. 1-28) 
where 
A = (2) 
mk =V, Hy (Via) Pat m #k 
(2) Zz 2 1 
ew Ip eo geal 
Mi | (3),, (2) 
V o ie (Vv, a) 
Bin = JOE g—- a7) eee 
1H Ow, a 
k=0 Vel o \Y, 
(Eq. 1-29) 
and 
N 2 
P = DOAN TL eed 
me ey (1:2) 
ay a (Eq. 1-30} 


The current as a function of distance along the antenna 


is then computed by 


ks 


| (2)=27 a H g(a,z ) (Om 201 


N 1 
=27n a > awa va) cos a a} 
C= 








M , 3 2 
Vo Lee a) mI Z 
=0 pel eallaoulyegdy) 


{Eq. 1-31} 
The input impedance is then given as the ratio of the gap 
voltage Vo to the current at the end of the gap 


Vo 


Zea 
"1 @ =d) {Eq. 1-32) 


B. SOFTWARE OVERVIEW 


The computer program was written in FORTRAN 77 and was 
designed to run on an MS-DOS compatible system with a numeric 
coprocessor and 460 kilobytes of available memory after the 
disk operating system in loaded. If the numeric coprocessor 
is not available, the source code could be recompiled to 
emulate the specific coprocessor calls. The user would notice 
a dramatic increase in run time, but as all calculations are 
done in single precision, the accuracy should not be affected. 

The ees code is divided into two separate programs. 


The first program, called MONO.FOR, is the main program that 


IE, 


calls a number of subroutines in SUBS.FOR. Variable names 
used in the program code closely follow those used in the 
theoretical development. Each variable is explicitly declared 
in each subroutine and array sizes are kept to a minimum to 
reduce required memory at run time. 

The first section of code in MONO.FOR is used to input 
the physical parameters using either an unformatted ASCII data 
file as shown in Appendix E or by entering the values from the 
keyboard. In each case, measurements may be entered in either 
wavelength or meters. If a data file is used, the first line 
must be either a "w" or "m", and may be either upper or lower 
case, to distinguish between wavelengths or meters. Unite 
number zero was used for all READ statements, which allows use 
of the program in a batch mode where all inputs and screen 
directed outputs can be redirected by standard DOS methods. 
However, the use of compilers other than MicroSoft Fortran 
(Versions 4.01 or 4.1) may produce compiler errors when using 
this designation. The user should consult the specific 
compiler manual for details on unit designations. 

Before prompting the user to identify the output file 
name, the program checks for two conditions that will produce 
inaccurate results. The first condition is when the upper 
plate height is a multiple of one half wavelength. Physical- 
ly, this condition will allow resonant cavity modes between 


the plates, introducing instabilities in the numerical 
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so luc ron ihe Mathematical "basis for this can be easily 
shown. When the upper plate height is a multiple of half the 
wavelength, v, will become zero for certain values of n. As 
can be seen in Table I, the value of a Hankel function ap- 
proaches -jo as the argument approaches zero. This test 
occurs in the main program and supplies earlier error detec- 
tion. The singularity is not checked for in the subroutine 
that computes the Hankel functions (HAN1) found in SUBS.FOR. 

The second condition that is tested is when the distance 
of the top ground plane above the antenna, designated qg, is an 
integer multiple of the antenna height. This is similar to 
the previous case in that a resonant area is present between 
the antenna and the upper ground plane. However, in this 
Gaccmmeene walue of 4, which appears only in the P,, matrix, 
becomes zero for multiples of the ratio hf to g. Figure 3 
illustrates the behavior of the Bessel function for a range 
of arguments [Ref. 4]. The zero order Bessel function, Jo, 
is finite for small arguments; however the first order Bessel 
function, J,, (found in the denominator of Equation 1-30) 
rapidly approaches zero as the argument approaches zero, 
Causing the errors. 

In an effort to reduce computation time at the expense of 
memory requirements, several arrays are dynamically loaded so 
that the values may be easily extracted to fill the system 


matrix A,, and the driving vector B,. These include the 


Poon 





Figure 3 Bessel Functions of the First Kind 


values for the Fourier moment integrations, the Bessel 
function values and the P,, matrix. The T matrix has dimen- 
sions of M, x N, and consists of terms from Equation 1-23. 
The I matrix has dimensions of M, x 3 where the column 
positions correspond to the three equations in Equation 1-18. 
The Ja and P,, matrices are loaded in the subroutine ALOAD and 
passed back to the main program. The Ja matrix is M, x 4 in 


size where the column positions represent the following 


values: 
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Column Position Value 


1 Un * Jo (Up a) 

2 Jy (uy, a) 

3 v,, * H,(v, a) 
4 (V,)° * Hy (Yq a) 


For large values of N,, N, and M,, about 60% of the total 
computation time is involved in filling the last temporary 
array, the P,, matrix. This array is dimensioned to N, x éM, 
and uses values from all three arrays for each location. By 
filling these arrays first, the computation time was decreased 
by an average of 85%, however the memory required increased 
by 300%! 

Once the system matrix and driving vector are stored in 
memory, it is a Simple task to solve for the unknown a,'s 
utilizing the back substitution method. First the matrix is 
upper triangulated by performing an L-U decomposition, with 
pivoting, on the A,, matrix using the FACTOR subroutine. Once 
this is accomplished, the SOLVE subroutine calculates the 
eigenvalues by back substitution. 

To maintain the convention that an outwardly travelling 
wave has a negative imaginary exponential argument, the CSR 
function was included in SUBS.FOR. Since not all compilers 
would return a negative imaginary value from the square root 


algorithm when the argument is a negative real value, the CSR 
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function insures the -j component of the square root of a 
negative argument is returned. 

The subroutines that compute the Bessel and Hankel 
functions for complex arguments and the matrix solution 
algorithms were supplied by Professor Michael A. Morgan of the 
Naval Postgraduate School. The Bessel and Hankel subroutines 
utilize the direct power series method for small values (less 
than or equal to five for the Bessel subroutine and six for 
the Hankel subroutine) and Hankel's asymptotic formula for 
larger values. 

The output is directed to a user specified file on any 
drive or subdirectory as long as the entire string is less 
than 25 characters. Longer strings will truncate the extra 
characters without warning. Printed output includes all input 
variables including the number of coefficients chosen for 
Regions I and II and for the gap voltage expansion. Length 
measurements are reported in both meters and fractions of a 
wavelength. The system eigenvalues are then printed followed 
by the position number, distance from the end of the gap and 
magnitude and phase of the computed currents, up to the upper 
plate height. Finally the input resistance and reactance is 
computed and printed to the file. Sample input and output 


files are listed in Appendix E. 
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Hl. VALIDATION 


In order to validate any computational algorithm, one 
must establish the accuracy obtained by comparing the computed 
results against either experimental data or another widely 
accepted software package that has been proven to be accurate. 
For this case, the logical choice was to use the Numerical 
Electromagnetics Code to calculate the input impedance and 
current distribution for a range of different physical 
structures and compare these results against those obtained 
by MONO. From this point, the program sensitivity will be 
examined for various ground plate spacings, antenna lengths 
aimee eradii. Chapter 4 will investigate the effect that 
reducing the number of modal expansion terms has on this 


accuracy. 


A. THE NUMERICAL ELECTROMAGNETICS CODE 


The Numerical Electromagnetics Code (NEC) is widely 
accepted as an accurate method of computing the antenna 
currents, input impedance and radiation patterns of antennas 
that can be modeled using wire or patch structures. NEC 
calculates these values by solving integral equations for the 
currents imposed on the structure by a voltage source or an 


incident electromagnetic wave. The structure is best modeled 
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by straight wire segments that should be less than ten percent 
of the operating wavelength and conforms to the physical 
geometry using a best fit approach. [Ref. 6, p. 3] 

The user supplies NEC with a formatted data file that 
specifies the location of the structure in the rectangular 
coordinate system, the number of segments used to compute the 
results, the driving source location and output parameters and 
the existence and characteristics of a ground plane in the XY 
plane. Because of the infinite number of images produced by 
a pair of perfect ground planes, the user must also construct 
enough images with the proper driving voltages that will yield 
consistent results. NEC then calculates the segment length, 
and uses this value as the gap distance between the ground 
plane and the base of the antenna. The currents are calcu- 
lated at the midpoint of each segment vice at the ends as is 
done in MONO. 

NEC uses two approximations in the calculation of the 
electric field, the thin-wire kernel and the extended thin- 
wire kernel. The thin-wire kernel can be used when the 
segment length to radius ratio is greater than 8, while the 
extended wire kernel is accurate when this ratio is at least 
2. Since NEC sets the gap distance equal to the segment 
length and accuracy is improved if the segment lengths 
(especially across a boundary) are of equal size, setting the 


gap distance equal to the antenna diameter forces the segment 
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length-to-radius ratio to be equal to 2. Additionally, this 
gap distance is consistent with the requirement for a constant 
electric field potential inside the gap region. Therefore, 
the extended wire kernel, which uses a uniform surface current 
along the segment length and assumes no variation of the 
currents along the @ direction, will be used for all valida- 
tion runs. These approximations also support the basic 
assumption of the TM mode dominating the currents on the 
monopole. For the conditions listed above, NEC's accuracy is 


within 1% of experimental results. [Ref. 6] 


B. COMPARISONS MADE 

The first goal will be to determine the minimum number of 
reflections that must be assembled in the NEC data set to 
obtain consistent results for a quarter wave antenna struc- 
ture. From this data, an appropriate number of reflections 
will be used that represents consistent results while reducing 
the calculation time for all future runs. This will be 
followed by demonstrating consistent results for various 
ground plane spacings. The upper ground plane was introduced 
to establish a closed region. The upper plate should, if far 
enough away from the antenna, have little effect on the input 
impedance and antenna currents for distances that are not 
exact multiples of half the wavelength. Again a value of / 


will be selected that represents reliable results but reduces 
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the calculation time required for MONO. The antenna radius 
will then be adjusted over the range of 0.005’ to 0.02, 
followed by variations in the antenna height for a selected 
radius. 

For all MONO calculations, the number of modes used in 
Regions I and II (N,, N,) will always be set to 60, while the 
number of modal expansion terms used for the gap field (M,) 
will be set to two times the upper plate spacing to gap 
distance ratio. Chapter 4 will further investigate the 
accuracy of the results as a function of the number of modes 
used. A preliminary check for consistent results can be 
obtained by evaluating the differences in the input impedance 
results. But these values may differ slightly since NEC 
calculates the current at the midpoint of each segment instead 
of at the ends of the segment as is done in MONO. Addition- 
ally, since the first segment used in NEC is the gap distance, 
the input impedance is actually calculated for a point not on 
the monopole but in the center of the feed region. In 
comparison, MONO calculates the input impedance at the base 
of the antenna. The effect this has on the comparisons will 
be discussed later in this chapter. Additionally, these 
differences may be emphasized since the magnitude of the 
current at the end of the gap region can be quite small so 
that even small differences in the calculated values will 


yield relatively high differences. Therefore, to better 
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assess the accuracy of the results, the calculated current 
distribution along the antenna will be compared graphically 


with that obtained from NEC. 


C. RESULTS 


Since the ground planes used to divide the area into 
regions create an infinite number of reflections, one must 
first examine the minimum number of reflections required to 
obtain convergence to the half space case, having no upper 
ground plane. For this set of runs, the antenna is 0.24) high 
With radius of 0.01. and the upper ground plate was placed 1.4 
wavelengths from the lower plate. Appendix F has a complete 
listing of some of the NEC data sets used for the validations. 
It should be noted that the distance between the lower plate 
and the center of the reflected dipole is equal to twice the 
plate spacing. This is due to the complete imaging, including 
the lower ground plane, of the physical structure from the 
upper ground plane. One can see from Figure 4 that adding a 
second reflection does not change the current distribution by 
any appreciable amount, therefore one reflection will be used 
for all further validations. 

Now that a baseline data set has been decided upon, the 
effects of various plate spacing, antenna heights and radius 
will be examined. The first item to look at will be the 


effect the plate spacing has on the current distribution and 
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the calculated input impedance. Table II shows the numerical 
results for an antenna of height 0.24 and radius of 0.01. 


placed between ground planes that are spaced 3.1, 1.4 and 0.81 


Table II RESULTS FOR VARIOUS UPPER PLATE SPACINGS 


a, MONO % Diff. 
a, Rin Xin R X 


44.63 13:535 42.65 13.14 4.64 1.60 
44.90 13.97 43.14 3 2 4.08 0.14 
43.51 14.84 44.67 13.84 2.66 6.74 





wavelengths. Figures 5 through 7 show the current distribu- 
tion for each of the above cases. While none of the results 
exceed acceptable limits, the case of /=1.4. is closer to the 
results computed by NEC and will be used as the baseline for 
further comparisons. It should be noted here that while the 
example values selected for the plate spacing have little 
effect on the result, the plate spacing can not be an integer 
multiple of a half wavelength. Likewise, the distance between 
the top of the antenna structure and the upper plate can not 
be a multiple of a half wavelength. Physically, these 
configurations would present the ability to support resonant 
modes that do not require a driving source. In reality, this 
condition will not exist as there is always a source of 
resistance which causes all modes to decay in time. However, 
the first condition will exist numerically where H,”’ ap- 


proaches negative infinity when vy, (Equation 1-11) equals 
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zero. By replacing the wave number k, with 27/\, then v, is 
zero when the plate spacing / is an integer multiple of a half 
wavelength. A similar condition will also occur when the 
distance between the top of the structure and the upper plate 
is a multiple of a half wavelength. In this case, J, found in 
the denominator of P,, (Equation 1-30) is zero when u, is zero. 
The program code will alert the user to both conditions and 
halt execution. 

The next variable to consider is the antenna height. For 
extremely short antenna structures, the current decays almost 
linearly along the wire from its initial value to zero. When 
the antenna is about one quarter wavelength long, the current 
distribution takes on the familiar cosine shape. As the 
antenna length approaches one half wavelength, the distribu- 
tion approaches a sine wave shape and the input impedance 
approaches infinity as the current approaches zero. 

Table III lists the results when the structure is placed 
between the two ground planes spaced 1.4 apart for various 
antenna heights and radii. In the case where the antenna 
height approaches a half wavelength, the percent difference 
in the input impedance values appears to be large. This can 
be attributed to the differences in the location on the 
structure where the input impedance is actually calculated. 
This disparity combined with the low magnitudes in general can 


lead to alarming discrepancies with the NEC results. 


3D 


Table III NEC, MONO RESULT COMPARISONS FOR VARIOUS hf AND a 
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Therefore, one must look at the current distribution to 
accurately determine the differences in the computed results. 
Figures 8 through 11 show that the differences in the 
magnitude of the current distribution are within two percent. 
Additionally, one could also compare the phase difference 
between the two methods to further establish an acceptable 
method of judging the accuracy. Figure 12 shows that the 
difference between the phase of the current as calculated by 
MONO and NEC is well within two percent along the entire 
length of the antenna. 

Another point that can be observed is the effect that a 
change in the antenna radius has on the results. Once again 
a review of the current distribution is required to accurately 
assess the effect any errors would have on the far field 


pattern; however one detail is evident. As the antenna 
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becomes thinner the difference between the two results is 
reduced. This effect can be attributed to the increased 
resolution in the NEC results as the radii and, therefore, the 
gap distance are reduced. Additionally, the approximation to 
a TM current distribution along the surface of the structure 
becomes more accurate with a thinner structure. 

In review, it has been shown that due to the differences 
in the methods, the best method to assess the accuracy of MONO 
using NEC is obtained by comparing both the current magnitude 
and phase distribution along the antenna surface vice the 
input impedance values. Using this criterion, MONO is 
consistently within two percent of the results obtained with 
NEC for a variety of structures. Chapter 4 will investigate 
the difference caused by reducing the number of expansion 
terms and attempt to identify a minimum number of terms 
required to obtain results that are consistent with those 


obtained here. 
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IV. CONVERGENCE CONSIDERATIONS 


In any realizable computing machine where there is a 
finite quantity of storage, the first question that must be 


answered is the number of terms required to obtain accurate 


results. If one uses too few terms, the results are er- 
roneous. However, using too many terms can often lead to 
inaccuracies as well. If the additional terms are small 


enough to be approximately on the same order as the precision 
of the computer, then the additional terms may actually 
increase the error. To compound the problem, the additional 
terms require additional memory and increase the time required 
to compute the results. In investigating this question, we 
will first examine the number of terms required for Regions 
I and II (N, and N,) followed by an analysis of the number of 
terms required to represent the electric field in the gap 


region (M,). 


A. MODAL TRUNCATIONS FOR REGIONS | AND Il 

In the previous chapter, the number of modes chosen for 
the two regions was the maximum allowed in the program, namely 
60. But is that enough? One could say that since the results 
compare favorably with NEC, the solution has converged. A 


more accurate approach is to increase the number of modes and 
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compare the results. However, due to the memory limit of 640 
kilobytes of addressable memory imposed by the compiler, the 
source code was recompiled on a 80386-specific FORTRAN 
compiler that allows for larger arrays and is only limited by 
the available memory of the computer being used. Figure 13 
compares the current distribution of a quarter wave antenna 
that has the upper plate height set at a large value (3.1)) 
for 60 modes and 100 modes. This additional information 
proves that the system has converged to an accurate solution 
with 60 modes and any additional modes are not required. The 
effect reducing this number has on the accuracy will be 


discussed later in this chapter. 


B. MODAL NEEDS FOR THE GAP REGION 
As mentioned earlier, the original approach was to expand 
the electric field along the surface of the antenna in one 


expansion between the ground planes as shown below 


V 
eens for O<zsd 
d 
Q, for d<zsh 
E, (a,2) = 
N, ] 
oe uC, Jn(tya) cos], forhszsl 
y n=0 : 


{Eq. 4-1} 
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However, this approach failed to yield consistent results and 
a review of the source code showed no errors in the logical 
design. Upon an examination of the electric field in Region 
I, it was discovered that a very high value for N, would be 
required in order to obtain convergence to the known field in 
the narrow feed gap. It was believed that significant errors 
in the electric field would cause inaccuracies in the result. 
The errors in the results were thought to be caused by the 
lack of terms used to describe the electric field. The 
electric field in the gap region can be thought of as a step 
function of short duration. In the frequency domain, func- 
tions of this type occupy a large spectral bandwidth, requir- 
ing many terms for the description. The large expected value 
of N, would necessitate the inversion of a very large matrix, 
thus making this approach inefficient. To avoid this problem, 
a separate expansion for the constant electric field was added 
to Equation 4-1 to form Equation 1-14. 

Once again the number of terms required was the main 
question. However, this field is known prior to the computa- 
tion and the number of terms can be estimated using standard 
Fourier analysis. MThe Fourier transform of a step function 
is known to be a sinc function with zero crossings at multi- 
ples of 1/d, where d is the duration of the pulse. Each lobe 
consists of impulse functions at intervals of the fundamental 


frequency, f,, which equals 1/T,, where T, is the period of the 
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pulse train. In this case, the upper ground plane causes a 
reflection of the original structure, but with opposite 
polarity, at a distance of 2/ from the lower ground plane. 
Therefore the number of impulses in each lobe of the sinc 
function equals 2//d. This can be investigated graphically as 
well. Figures 14 to 19 (pages 48 to 50) illustrate the field 
computed for a structure with a gap distance (d) of 0.06 
meters and upper plate height of 9.54 meters and fed by a 
constant one volt source for various numbers of terms in the 
series. As can be expected, when the number of terms are not 
an integral multiple of //d, excess oscillations occur. It is 
also important to note that as the number of terms increase, 
the calculated curve approaches the ideal value at the 
discontinuities, while the pointwise error along the rest of 
the curve increases. However, the error in the least square 
sense decreases. This is known as the Gibbs phenomenon. Near 
a discontinuity for a large number of terms, the overshoot is 
approximately 8.95 percent of the desired value, or -18.17 
volts [Ref. 5 p. 557]. This value is reached when the number 
of terms equals 3//d as seen in Figure 19, therefore addition- 
al terms will not reduce the error at the discontinuities. 
From these graphs, one can see that the number of terms used 
to represent the electric field in the gap region should be 
an integer multiple, less than or equal to three, of the ratio 


of upper plate height to gap distance. An upper plate height 
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to gap distance ratio of 250 was considered to be the maximum 
desired. Therefore to insure the code would fit easily into 
the 640 kilobyte limit, the number of terms used for this 
expansion was set at 2//d. This was the value used for all 


the validation calculations completed in Chapter 3. 


C. EFFECT OF REDUCING N, AND N, 

If the number of terms can be reduced and still provide 
consistent results, one can then decrease the size of the 
internal arrays and reduce the computation time. For example, 
by reducing N, by 50%, the computation time is reduced by 
almost 48%, and when both N, and N, are reduced by 50%, the 
computation time is reduced by 71.5%! However, the current 
distribution along the antenna is unknown, making it difficult 
to calculate the number of coefficients required for consist- 
ent results based on the input parameters. The results 
obtained by independently varying N, and N, are shown in Table 
IV and are compared to the results obtained in Chapter 3 for 


an antenna with the following physical parameters: 


Antenna Height (h) 0.24) 
Antenna Radius (a) 0.01) 
Gap Distance (d) 0.02. 
Plate Spacing (/) Zak TX 
Number of terms for 217 


Gap Region (M,) 


Sl 


Table IV COMPARISON OF RESULTS AS N, AND N; VARY 





As expected, the difference in the input impedance increases 
as the number of terms is reduced. However, the interesting 
point is the effect that the apparent relationship between N, 
and N, has on the result. If N, and N, are equal, acceptable 
results are obtained for as low as 30 terms, while one must 
use at least 50 terms if the values are not equal, to get 
comparable results. A mathematical explanation for this 
relationship is not yet understood and will be investigated 
in future thesis efforts. Figure 20 represents the current 
distribution as N, and N, vary independently compared with the 


result when N, and N, equal 60. 
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Figure 20 Current Distribution as N, and N, Vary 


D. EFFECT ON ACCURACY AS M, VARIES 

As seen in the previous section, one can reduce the 
memory requirements and computation time by reducing the 
number of coefficients used in Regions I and II. Using 
Fourier analysis, the user can calculate a value for the 
number of terms required to accurately represent the electric 
field along the gap (M,). Since this value is much larger 


than N, and N,, one could dramatically reduce the required 
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storage and speed up the computations by reducing M, to the 
same order as N, and N,. 

To demonstrate the sensitivity of the input impedance to 
the number of terms used for the electric field along the gap, 
one can compare the results against those found when N, and N, 
are 60 terms and M, is computed as 2//d. Table V lists the 
results for same antenna structure described on page 51. As 
expected, the errors in the calculated input impedance begin 
to increase as the number of terms are reduced. However, this 
difference is not linear for all values of N, and N,. In fact, 
Figure 21 suggests an optimal number of M, terms, for given 
N, and N,, exists which results in smaller differences between 


the two cases. Since M, is a function of the gap distance, 


Table V MONO SENSITIVITY TO (M,) FOR GIVEN N, AND N, 
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the input impedances were recalculated when the gap distance 
is reduced to 0.01\ for selected values of N, and N, as shown 
in Figure 22. While this clearly demonstrates that an optimal 
number of terms may exist that produces results which approach 
those for large values of M,, N, and N,, a clear relationship 
between the physical structure and these values are not yet 
fully understood. 

But how can the number of terms be reduced by almost 91 
percent, and yet the difference is no greater than 5.13 
percent? First, the basis for selecting the number of terms 
used was based on obtaining accurate representation for the 
electric field along the antenna. However, since the far 
field patterns are derived directly from the currents, it is 
the current that must remain consistent and not necessarily 
the electric field along the surface of the antenna. From 
Equation 1-31, we see that the current is derived from the H,, 
or the first derivative of the potential field with respect 
to p. However, the electric field along the z direction is a 
function of the second derivative of the potential with 
respect to z. Since the derivative process increases the 
noise, or errors, it is evident that small differences in the 
potential field may yield wide variations in the electric 
field while the magnetic field, and therefore the current 
distribution, may remain within accepted accuracy limits. To 


prove this point, Equation 4-2 shows an expression that 
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approximates the electric field along the surface of the 


antenna by taking the second derivative of Equation 1-31 with 


| 


respect to z. 
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(Eq. 4-2) 
This equation can now be used to visualize the approximate 


change in the electric field along the antenna for a given 
number of modes M,. As shown in Figure 23, the change in the 
magnitude of the current distribution when M, is changed from 
217 (which equals 2//d) to 30 is negligible. However, Figure 
24 shows the curves obtained from Equation 4-2 for the same 
parameters as used in Figure 23. As one can see, the field 
is accurately represented near the end of the wire, however 
near the gap region, there are an insufficient number of 
coefficients to accurately compute the fields in this region. 
This refutes the original assumption that the electric fields 
must be accurately represented to guarantee an accurate 
current solution. Additionally, since M, can be on the same 
order as the values used for N, and N,, then the coefficients 
used to describe the gap field and the field along the 


structure can be combined into one coefficient. The next 
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chapter will explore the theory and the differences this new 
approach has over calculating a separate set of expansion 


coefficients for the gap region. 
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Figure 23 Current Distribution as Ml Varies (h=0.24)) 
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Figure 24 Approximated Electric Field Distribution 
(Calculated using Equation 4-2) 
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V. A SIMPLER APPROACH 


During the early stages of the software development, it 
was decided that a separate expansion for the constant 
electric field in the gap region was required to obtain 
consistent results for the input impedance and current 
distribution along the structure. This decision was based on 
the belief that large errors in the total electric field along 
the structure would result in large errors in the current 
distribution. However since the current distribution is 
derived directly from the magnetic field, and is relatively 
insensitive to the electric field, the number of terms used 
for the gap region expansion was reduced to the same order as 
those used for N, and N, without introducing appreciable 
errors. Noting this, the original method of describing the 
system was reevaluated such that one could reduce the required 
memory and speed up the calculations even further by combining 
the a, and b, coefficients used in Equation 1-10 (page 11) into 
one set of terms and by deleting the separate expansion along 
the gap region entirely. This chapter will work backwards 
from the technique presented in Chapter 2 and will note the 
changes in the system equations resulting from this approach. 
The resulting code, which will be referred to as MONO8, will 


then be compared to the results obtained for the cases with 
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the largest differences between NEC and the code using the 


separate expansion (called MONO7). 


A. MODIFYING THE TECHNIQUE 


As was discussed in Chapter 2, the previous approach used 
two expansions in Region I where p 2 a. The first expansion 
sets the tangential electric field to zero from the lower 
ground plane to the top of the antenna. A second expansion 
sets the tangential electric field to - V,/d between the lower 
ground plane and the feed point of the antenna and sets E, to 
zero from this point to the upper ground plane. Since the 
Fourier transform of a sum is the sum of the transforms [Ref. 


7), we can combine the expansions for the potential field in 


Region I as shown below 


1 


N | 
y (p,z)= » said, 0 (Vp cos z : 
nay (Eq. 5-1} 





where a, equals the sum of the a, and b, coefficients used in 
Chapter 2. Again we match the magnetic fields across the 


boundary between Region I and II such that 
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(Eq. 5-2) 
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By utilizing the orthogonality property of the Fourier moment 


integrations used in Equation 1-15, the system equation can 


be expressed as 


N, N 4 
di el eete 1) of tmaeer sa 2 2 (1) 
So vpapHy (ya) Deen nat Vm Hg (Vn @) Im 
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= me 


(Eq. 5-3} 


for m= 0 to N,. This can be expressed in matrix form as 


Ny 
> Amc = By for m =0,N, 
k=0 
(Eq. 5-4) 
where A,, is given in Equation 1-29. As the b, coefficients 
are now included in the a, terms, the driving vector, B,, is 
reduced to 


|b a 


m 


@Vo 
m 

t (Eq. 5-5} 
Once the unknown a, coefficients are found, the current along 


the antenna, as a function of z, is now given by 


NG 


I (z) = 27a SS va, H 3 (Vv, a) cos 2] 


n=0 
{Eq. 5-6) 

B. SOFTWARE MODIFICATIONS 
As seen in Equations 5-3 to 5-6 above, removing the 


separate expansion for the electric field along the gap region 
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requires changing the software code for only the driving 
vector B, and the expression for the current as a function of 
position along the antenna structure. The first advantage 
realized was the reduction in the amount of memory storage 
required to run the progran. Since it was shown in the 
previous chapter that the system had converged when the number 
of modes in Regions I and II were set at 60, independent of 
the number of modes in the gap region, it was not necessary 
to dimension any array greater than 60. However, we can now 
reduce the size of the I, T, Ja, and P,, matrices by almost an 
order of magnitude, such that the amount of memory required 
to run the program has decreased by 75 percent to less than 
104 kilobytes. Additionally, the execution time of the 
program has decreased by an average of 35 percent for the runs 
presented here. 

But does this method yield results that are comparable to 
the results when the separate expansion is used? We will 
examine the differences between MONO7 and MONO8 for various 
plate spacings, antenna heights, as well as the effect that 
independently varying N, and N, has on the results obtained 


earlier. 


C. RESULTS 
As a first test of the accuracy for MONO8, we will 


examine the differences in the calculated input impedance 
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between MONO7 and MONO8 as the ground plane spacing is varied. 
The three structures used in Chapter 3 for Table II were rerun 


where N, and N, were set to 60. Table VI indicates that the 


Table VI COMPARISON OF MONO7, MONO8 FOR VARIOUS ! 


13.14 43.09 13.12 1.03 
13.95 43.15 13.96 0.02 
13.84 44.68 13.82 0.02 





results obtained from the new method differ only slightly from 
the results obtained when the gap expansion is used. 

The next test of the accuracy considered different 
antenna heights. In this case, the comparison is made between 
the structures that resulted in the highest difference between 
MONO7 and NEC, or when the antenna radiuS was 0.02. Table 
VII shows the results of this comparison for a plate spacing 
of 1.4. Once again, the calculated input impedances differ 
very little from those obtained using the separate gap 


expansion. 


Table VII COMPARISON OF MONO7, MONO8 FOR SELECTED A 


219.44 |-121.34 





64 


As noted in Chapter 4, the results are more consistent 
with the solutions obtained from NEC when N, equals N,. Figure 
25 illustrates that this same phenomenon occurs for MONO8. 
Further investigation and study is required to understand the 
relationship that the number of modes in Regions I and II have 


on each other. Once this is understood, it may be possible 


Compar ison of the Current Distribution 


aon For selected N1,N2 respectively 


.028 
.026 
.024 
D022 
.020 
.0168 
.016 
.014 


012 


fan 
ify) 
@ 
WS 
® 
Qa 
= 
© 
Ww 
~ 
(= 
v 
= 
= 
= 
© 


.010 
.008 
.006 
. 004 
. 002 


.000 
U2 a4 


Distance along Antenna Cmeters) 
+ 20,30 ° 30,20 


Figure 25 Current Distribution using MONO8 as N, and N, Vary 





to derive an expression for the minimum number of modes 
required to obtain results within a specified accuracy range. 
This expression could then be used to allow the software to 
dynamically adjust the number of modes required to obtain 


accurate results at run time. 


65 


In all the above cases, the results were compared to the 
solution using MONO7 where M, was chosen as 2/l/d. However, 
Figures 21 and 22 (page 56) indicates that there is an optimum 
number of modes for the gap region to reduce the differences 
between the solutions of MONO and NEC. Table VIII shows that 
MONO8 returns the same result for the input impedance as MONO7 
did where M,, N, and N, were equal. While this result may not 


be the closest to NEC's results, the differences are consis- 


Table VIII COMPARISON OF MONO8 AND MONO7 AS M, VARIES 


MONO8 MONO7 R/X 
N1=N2 R/X M1=60 M1=50 M1=40 M1=30 
43.31 43.31 43.02 43.39 42.44 
13.87 13.87 13.85 73.95 14.14 
42.81 42.59 43.12 42.13 
13.44 13.42 13.53 13.72 
30 40.70 41.61 41.66 41.72 40.70 
12.64 12.32 12.32 23.39 12.64 


tently less than ten percent. 


















In conclusion, consistent results are obtained without 
the separate expansion for the gap region. While these 
results may not offer the closest possible solution to that 
available by NEC, the small difference in the solution 
obtained without the extra terms is more than offset by the 


significant reduction in required memory and calculation time. 
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When the relationship between the physical structure and the 
Minimum number of modes that are required to obtain the best 
possible solution is understood, then the software can be 
modified so as to select the number of terms required at run 
time. Until then, MONO8 should be the method of choice as 
its results are more predictable and is capable of returning 


the result within acceptable accuracy limits without "trial 


and error". 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


We have shown that one can calculate the current dis- 
tribution along a monopole using multiregional cylindrical 
harmonic expansions that are consistent with the results 
obtained from the Numerical Electromagnetics Code (NEC). 
While this software package is not intended as a replacement 
for NEC, it is the beginning step in the development of 
software directed towards solid antenna structures that are 
rotated about a symmetrical axis with a simple and flexible 
input structure. The particular software that has been 
developed has the following characteristics: 

- Simple input data structure. All inputs entered direct- 
ly, via unformatted ASCII file, or redirected using 
standard DOS redirection codes allowing batch processing 
of multiple runs. 

¢ Very flexible for a given monopole antenna structure. 

* Can be readily modified to allow for a homogeneous 
dielectric coating on the surface of the antenna. This 
coating could be used to modify the surface currents and 
the far field radiation pattern. 

As currently written for the simple monopole, the user is 
able to modify all the physical parameters of the structure 


including driving source frequency. With some modifications, 


including the incorporation of a finite element algorithm in 
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one sub-region, the user may also include the effect of a 
nonhomogeneous dielectric coating so that the far field 
radiation pattern may be tailored as desired. 

One difference that should be noted between the techni- 
ques used by NEC and MONO is the location used to calculate 
the input impedance. Due to the NEC's segmentation of the 
entire structure, the input impedance is actually calculated 
in the center of the gap region. However, MONO calculates the 
input impedance based on the value of the current at the base 
of the antenna. This difference may lead to large discrepan- 
cies in the calculated input impedances if the value of the 
current is changing rapidly at the end of the gap region, as 
was observed for monopoles whose height approaches a half 
wavelength. Other differences include: 


¢ MONO requires longer computation time for larger struc- 
tures than does NEC. 


« MONO Requires larger available memory. 


¢ MONO currently has limited error detection and correction 
within the subroutines. 


¢ The current inability to calculate a number of internal 
parameters requires a "trial and error" approach to 
determining the minimum number of modes required for 
proper convergence. 


- User is required to make extensive modifications to 
MONO's code for various types of antenna structures. 


- MONO is a highly specific, concept validation platform 
vice a generally applicable software package. 
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B. RECOMMENDATIONS 

Continued research is required in order to understand the 
effect that the apparent relationship between N, and N, has on 
the accuracy of the result. Additionally, there appears to 
be an optimal number of terms required for the gap region that 
results in an apparent improvement of the result when compared 
to those obtained using NEC. These unanswered questions point 
out the need for additional analysis of the relationship 
between the number of terms required to represent the electric 
field in a specific region and the overall accuracy of the 
final result. When this is accomplished, modifications should 
be made to allow the software to calculate the optimal number 
of terms required for each expansion that would yield results 
in the shortest time possible for a specified accuracy range. 

Another modification could be made to the code that would 
reduce the amount of memory required to run the program. This 
can be accomplished by temporarily storing only one row of the 
P,, matrix at a time, vice storing the entire matrix. This 
approach could reduce the amount of required memory by as much 
as 50 percent. 

To better assess the benefits of this method over the one 
used in other computational schemes, experimental data should 
be compared to the computed results for a variety of struc- 


tures. In this way, a basic understanding of when this method 
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best approximates the actual performance characteristics of 
an operational antenna can be obtained. Once these areas are 
better understood, a collection of routines for various 
physical structures should be developed using the basic 
technique developed in Chapter 2. This collection could 
include top loaded monopoles, and conical and other structures 


that feature symmetry about a axis of rotation. 
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APPENDIX A 


SOFTWARE FLOW CHART FOR MONO WITH GAP FIELD EXPANSION 
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Open 
Filename 


Write tnput 


Varlables to File 


Load Tran 
Matrix 





write 
O tru 
a Ee eS a SE = 
Solve for 
pa F ae ert 
Ja Matrix unction of Position 
ALOAD 
Subroutine 


Load Pmk 


Write Results 
Matrix 


to 
Output File 


Load System 
Matrix (Amk) 


Compute 
Input 
impedance 





Load Driving 
Vector (Bm) 


Write Result 
Cal) 
Factor 
Subr out ine 


Call 
Solve 


to 
Output File 








Subroutine 
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ONONND ANA YANAANAAAIANAIAANANANANAAAANAAAAAANH 


OooIaaqnrianrnan 


APPENDIX B 


SOURCE CODE FOR MONO WITH GAP FIELD EXPANSION 


PROGRAM MONO 
(Revision 7.3) 


Computing Currents on a Monopole Between Parallel 
Plates using Cylindrical Harmonic Expansions in 2 Regions. 


Links with SU8S4.0BJ. Created 05/87 LT. R.C. Hurley 
ae Variables© <---> -c----- 

FMHZ = Frequency in MHz 

h = Monopole height in Meters 

3 = Monopole radius in Meters 

d = Distance between lower ground plane and feed point 
in Meters 

l = Distance between ground planes in Meters 

q = height from top of monopole to upper ground plane 
in Meters 

J = Imaginary One 

pi = Value of pi (3.1415927) 

KO = 2*pi*F/c 

where F = frequency of operation in Hz 
c = Speed of Light 

K2 = k0**2 

N1,N2 = Number of Coefficients in each of the two regions 

NBIG = Size of the arrays used in FACTOR and SOLVE 

Unit 6 = Screen Output 

Unit 7 = Main Output File to Disk (Name selected by User.) 


xeeek Variable Definitions ***** 


INTEGER k,m,M1,Mmax,mon,n,ns,N1,N2,Nbig,NI ,Nmax,P(61) 

REAL a,ai,Am,An,Bl,c,CI,CM,CP,CR,d,di,DeltaZ,Denom,eps, FMHZ,h,hi 
REAL 1(0:500,1:3),KO,K2, |, li, Lambda, Pi,g,Rin,S 

REAL T(0:500,0:60),w,Xin,y, yhat 

COMPLEX Amk(61,61),8(61),C1,C2,1Z(61), j, J0,J1, Ja(0:500, 1:4) 
COMPLEX Pmk(0:60,0:500),Sum,Sum1,Y0,Y1,Zin 

CHARACTER*25 Fname 

CHARACTER*1 Bell 

CHARACTER*1 HOW 

CHARACTER*1 INPUT 

CHARACTER*12 INname 

COMMON/BLk1/1,j,N1,Pi,d,Nbig 

COMMON/BLk2/K2,\,Mmax,N2,q,a,T 


RARAEKEKKERAERAREEKREKRERRERKREEREKKEE 


. Main Program Code I 


REKKKEEAKEEKAKEEREEEEEREREKRKKKEKRKRKRKK 
RERAEEKARKEKEKKKEKKRERRKEER 


i Input Section = 
RAREREEKEREEKKKKKKKKKKRKKE 


Bel L=CHAR( 7) 
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Nbig=61 

c=3.0E+08 

WRITE(6,*) ‘Do you wish to use an input data file [Y/N]?! 
READ(O,100) HOW 

IF CHOW.EQ.’y’ .OR.HOW.EQ.’Y’) THEN 


+ ‘m’) 


WRITE(6,*) ‘Enter input data file name with extension :/ 
READ(0,100) INname 
OPEN(Unit=4, File=INname) 
READ(4,100) INPUT 
IF (CINPUT .EQ. ‘w’) .OR. CINPUT .EQ. /W’) .OR. CINPUT .EQ. 
-OR. CINPUT .EQ. /M’)) THEN 
READ(4,*) FMHZ,hi,ai,di,li,N1,N2,M1,NI 
ELSE 
WRITE(*,*) Bell 
WRITE(6,*)’Error detected in first line of data file ’ 


+, INname 


WRITE(6,*) ‘First line must be either M or W for meter 


+ values or wavelength values.’ 


ELSE 


+/W)?! 


STOP 
ENDIF 


WRITE(6,*) ‘Enter data in either ([MJeters or [WJavelength [M 


READ(0,100) INPUT 

IF CINPUT.EQ.’m’ .OR.INPUT.EQ./M’) THEN 
WRITE(6,*) ‘Enter Frequency in MHz: / 
READ(O,*) FMHZ 
WRITE(6,*) ‘Enter Monopole Height, h, in meters: 
READ(O,*) hi 
WRITE(6,*) ‘Enter Monopole Radius, in meters: 
READ(O,*) ai 
WRITE(6,*) ‘Enter Feed Gap, d, in meters: 
READ(O,*) di 
WRITE(6,*) ‘Enter Upper Plane Height, | (meters): ‘ 
READ(O,*) Li 

ELSEIF (INPUT.EQ.’w’ .OR.INPUT.EQ./W’) THEN 
WRITE(6,*) ‘Enter Frequency in MHz: ° 
READ(O,*) FMHZ 
WRITE(6,*) ‘Enter Monopole Height, h, as factor of wav 


+elength : / 


READ(O,*) hi 
WRITE(6,*) ‘Enter Monopole Radius, as factor of wavele 


+ngth : / 


READ(0O,*) ai 
WRITE(6,*) ‘Enter Feed Gap, d, as factor of wavelength 


READ(0,*) di 
WRITE(6,*) ‘Enter Upper Plane Height, l, as factor of 


+wavelength : ’ 


READ(O,*) li 
ELSE 
WRITE(*,*) Bell 
WRITE(6,*) ‘You must enter either an M or W. 


+Reenter.. / 


GOTO 2 
ENOIF 
WRITEC6,*) ‘Enter # Region 1 Coeffs. (N1 .LE. 60): : 
READ(O,*) N1 
WRITE(6,*) ‘Enter # Region 2 Coeffs. (N2 .LE. 60): 4 


READ(0O,*) N2 
WRITE(6,*) ‘Enter # Gap Region Coeffs.(M1 .LE. 500): ’ 
READ(O,*) M1 
WRITE(6,*) ‘Enter No. Points for I(z) including ends: / 
READ(O,*) NI 
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AMAaAANAAnNAANN 


AOA AAAN 


NAAaAAAAaNAA 


ENDIF 


IF (M1.GT.500) Then 
WRITE(*,*) Bell 


ENDIF 


WRITE(6,*) ! 
WRITE(6,*) ! 
WRITE(6,*) ! 
STOP 


Ratio l/d too large. Adjust either: ’ 
a) l or d such that 2l/d < 500’ 
b) Array sizes of Im,T,Ja’ 


REKKKKEKKKKKKKKEKKKKKEKKKKKKKKKKKKKEKKKKKKKKKKKKEK 


* * 
- Calculate the Physical Parameters : 
* in both Wavelenght and Meters = 


* 


* 


RARKKKKKKKEKEKRKEKEKKRKEKAKKKKKEKKEKKEKEKEKKKKKKKKKKKKKEK 


Lambda=c/( FMHZ*10**(6)) 
IF CINPUT.EQ.‘w’ .OR.INPUT.EQ./W’) THEN 


ELSE 


ENDIF 


h=hi*Lambda 
a=ai*Lambda 
d=di*Lambda 
l=li*Lambda 


h=hi 
hi=h/Lambda 
a=al 
ai=a/Lambda 
d=di 
di=d/Lambda 
l=l1 
li=l/Lambda 


RRRRARKRAKKKKAKKKKEKEKKEKEKREKKEKKKEKKKKEKKKAKKKKKK 


® 
* 
® 
* 
Kakkkk&kkK 


* 


Is Upper Plate Height an * 


Even Multiple of Wavelength? sl 
*® 


RRR KKKKKKKKKEKKKKKKEKEKKKEKKEEKEKKEKKK 


IF (MOD((L*2),Lambda).LT. 1E-05) THEN 
WRITEC*,*) Bell 


WRITE(6,*) ! 
WRITE(6,*) ! 


f 


ERROR! Upper Plate Height 1s almost an even mul 


+tiple of half the wavelength. ’ 


WRITE(6,*) / 


Ensure l is NOT within 0.00001 of an integer mul 


+tiple of half the wavelength. ’ 


ENDIF 


STOP 


kkkkkkke 
* 
x 
* 
* 
kkkkkkkn 


RARKKAKKKKKEKKREKEKEKKKKERKKEKEKEKKEKKKKEKK 
* 


Is Upper Plate Height an id 
Even Multiple of Antenna Height? ii 
* 


RREKKKEKKKKEKKEKEKEKKEKEEEEKKEKREKKEKEKKKEKK 


E=(10*h)- (ANINT((10*h)/L)*L) 
IF ((MOO(L,h)).LT.1E-05) THEN 
WRITE(*,*) Bell 


WRITE(6,*) / 
WRITE(6,*) ! 


t 


ERROR! Upper Plate Height is almost an even mul 


+tiple of the antenna height.’ 


WRITE(6,*) / 


Ensure l is NOT within 0.00001 of an integer mul 


+tiple of h (in wavelength).’ 
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STOP 
ENDIF 
WRITE(6,*) ‘Enter output file name and extension: ‘ 
READ(0,100) Fname 
OPEN (Unit=7, File=Fname) 
WRITE(7,210) FMHZ 
WRITEC7,220) h,hi 
WRITEC7,230) a,ai 
WRITE(7,240) d,di 
WRITEC7,250) L,Li 
WRITEC7,*) ‘Region I Coeffs. (N1): ’,N1 
WRITEC7,*) ‘Region II Coeffs. (N2): ‘,N2 
WRITEC7,*) ‘No. Gap Coeffs. (M1): /,M1 
WRITEC7,*) ‘No. Points for I(z): ‘NI 
j=(0.0,1.0) 
Pi=3.1415927 
KO=Pi*FMHZ/150.0 
K2=K0O*KO 
g=l-h 
Nmax=MAX0(N1,N2) 
Mmax=MAX0(Nmax,M1) 
WRITE(6,*) ‘Loading the Moment Matrices ..... a 
WRITE(6,*) 4 / 


C kkkkkkkkkkkkkkkkkkkk&kkkk&kkkkkkkkkkkkkkkkkkkkkxk 
C * * 
C : Evaluating Moment Integrations td 
c * Over Region h < z < lL. . 
C * ey 
C kkkkkkk&kkkkkkkk&kkkkkk&kkkkkkkkkkkkkkkxkkkkkkkkkK 

T(0,0)=q 

DO 11 k=1,N2 

T(0,k)=0.0 


11 CONTINUE 
DO 22 m=1,Mmax 


Am=m*Pi/l 
DO 22 n=0,N2 
An=n*Pi/q 
IF(CABS(Am-An)).GT.1.0E-3) THEN 
Denom=Am* Am- An*An 
S=SIN(CAm*h) 
T¢(m,n)=-Am*S/Denom 
ELSE 
mon=m+n 
T¢m,n)=0.5*%q*((-1)**mpn) 
ENDIF 
22 CONTINUE 
WRITEC6,*) / Moment Matrix T Loaded. ’ 


WRITE(6,*) / ’ 


c KkkkKKKKKKKK&KKkkkkkkKKkKkKkkKKKKKKKKKKkkk&kKkkkkkxkk&kE 
C * & 
Cc ol Loading the Im Matrix for the : 
c . Largest Number of Coefficients . 
c * e 
€ kkekkkkkxkkkkkkkkkkkkikkkkkkkkkkkkkkxkxkkkkkkkkkkkkk: 


00 33 m=0,Mmax 
IF(m.EQ.0) Then 
I¢m, 1)=l 
I(m,2)=q 
I(m,3)=d 
ELSE 
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33 


agaAAaAaANn 


25 


I¢m,1)=l/2.0 
I(m,2)=q/2.0 
Am=m*Pi/| 
1¢m,3)=SINCAm*d)/Am 
ENDIF 
CONTINUE 
WRITE(6,*) ' Moment Matrix I Loaded.’ 
WRITE(6,*) ° 
WRITE(6,*) ‘Loading System Matrix .....’ 
CALL ALOAD( Ja, Pmk, Amk) 
WRITE(6,*) ‘System Matrix Amk Loaded.’ 
WRITE(6,*) ° ? 
WRITE(6,*) ‘Loading Driving Vector ..... f 
CALL BLOAD(FMHZ,M1,Pmk, Ja,B) 
WRITE(6,*) ‘Driving Vector B Loaded.’ 
WRITE(6,*) ¢ 
WRITE(6,*) ‘Solving the Linear System .. Please Wait .... ’ 
Ns=N1+1 
CALL FACTOR(Amk,P,Ns,Nbig) 
CALL SOLVE(Amk,B,P,NsS,Nbig) 
WRITE(7,*) ‘Listing the Expansion Coefficients’ 
WRITE(6,*) ° Listing the Expansion Coefficients ...../ 
WRITE(7,*) ‘N and B(N): ? 
DO 44 N=0,N1 
WRITE(7,*) N,BCN+1) 
CONTINUE 


RERKAEKEAKKEKEKEKEKEKKKEKKKEKKKEKEKKKEREKKKKEKEEK 


* * 
. Computing Current from the x 
. Expansion Coefficients = 
* x 


RERAAKAKKAKKAKKKEEREKEKKKEKEKEEKEKEKEKEKRKEKKKEKKKE 


WRITE(6,*) ' ? 

WRITE(6,*) ‘Computing current on the monopole..... ‘ 

WRITE(7,*) ‘Listing the Current on the Monopole’ 

DeltaZ=(h-d)/(NI-1) 

WRITE(7,*) ‘Position Distance 1(z) (A)! 

WRITE(7,*) ’ Number (Meters) Mag Phase’ 

eps=8.854E-12 

w=2 .0E6*Pi *FMHZ 

yhat=w*eps 

y=d 

k=1 

Sum=(0.,0.) 

DO 66 n=0,N1 
Sum=Sum+ (B(n+1)*Ja(n,3)*cos(n*Pi*y/l)) 

CONTINUE 

Sum1=(0.0,0.0) 

DO 77 m0,M1 
Sum1=sum1+((1(m,3)*Ja(m,3)*COS(m*pi*y/l))/(Jacm,4)*1(m, 1))) 

CONT I NUE 

1Z2(k)=2.0*P1*a* (Sum- j*yhat*Sum1/d) 

CM=CABS(12Z(k)) 

CR=REAL(IZ(k)) 

CI=AIMAG(IZ(k)) 

CP=180.*ATAN2(CI,CR)/pi 

WRITE(7,200) k,y,CM,CP 

y=y+Del taZ 

k=k+1 

IF(y.LE.(h+1.E-6)) GO TO 55 

WRITE(6,*) ‘Complete.’ 
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aqaanaan 


onqonInann 


100 


103 
200 


210 
220 
230 
240 
250 


RARAKAAKAARKAEREKAREREEEREEREEEEEEEEEEREEREEREEKEEEK 


* 


5 Computing Input Impedance 


* 


* 
* 
*& 


RHEE EEEEREREEEEREREEREEEEREEEREREEE 


Zin=(1.0,0.0)/12(¢1) 

RiN=REAL(Zin) 

X iIN=AIMAG( Zin) 

WRITE(7,*) ‘Input Resistance = ',Rin 
WRITE(7,*) ‘Input Reactance = ‘,Xin 


RAARAKEAKKAKEKEKEKKAKEKEKEREREKEREEEEEEEREEKEKEKKREE 


* 


* 1/0 Format Statements 


* 


* 
* 
* 


ARARAAAAEEEAAKAEEKREEEERAEKKEREEEEREEEEEREEKEKEKE 


FORMAT (A) 
FORMAT(15) 
FORMAT (E12.3) 


FORMAT (2X,13,5X,F8.3,5X, 1PE12.3,5X, OPF8.2, 1PE12.3,5X,OPF8.2, 


+1PE12.3,5X, OPF8.2) 
FORMAT(’ Driving Frequency, ',9X,F7. 
FORMAT(’ Monopole Height, 1 9X, F7. 
FORMAT(’ Monopole Radius, ‘,9X,F7. 
FORMAT(’ Feed Gap Distance, ',9X,F7. 
FORMAT(’ Upper Plane Height,’ ,9X,F7. 
STOP 
END 


5,! 
Ae 
3 / 
3,7 
bi 


MHz’) 

meters 
meters 
meters 
meters 


qe, 


(',F5.3,' Lambda)’) 
(*,F5.3,/ Lambda)’) 
(',F5.3,' Lambda)’) 
(1,F5.3,/ Lambda)’) 


MAAAANAANANAAANANANANANAAANANANANANANAANAANINANANANNANANANNANAANAAAMN 


ANAAAANANANANAANAAN 


Subroutine Source Code for MONO with Gap Field Expansion 


SUBROUTINE ALOAD(Ja,Pmk, Amk) 
(Revision 7.3) 


RHEAAEEKAKAKEKEKREEAKRERAAEAEKEEKEEEEEKEKEKAKKKS 


* * 
* Variable Definitions * 
* Sd 


kkekkkkkkkkkkkkkkkkkkkkkkkkhkkkakkkkkhkhkakkhkkh 


l Distance between ground planes 


Radius = Radius of monopole 

q = Distance between top of monopole and upper plane 
d = Antenna feed gap 

K2 = KO**2 (From main program) 

T(m,n) = Moment integration result as function of m and n 
I(m,x) = Moment integration result of order x given m 

j = Imaginary One 

Zz = Intermediate value for Bessel Subroutine 

un, vn = Eigenvalues 

JO = Bessel Function of the First Kind of Order 0 

Ji = Bessel Function of the First Kind of Order 1 

HO = Hankel Function of the Second Kind of Order 0 

H1 = Hankel Function of the Second Kind of Order 1 
SUM = Summation over n for the Amk matrix 

$1 = Individual term used in SUM 

m = Column position of matrix Amk 

k = Row position of matrix Amk 

Ni = Number of suddivisons in Region 1 

N2 = Number of subdivisons in Region 2 

Amk = Complex Matrix value that describes the system 


This SUBROUTINE loads the complex A matrix that describes the system 
of size Ni by Ni. Maximum size of matrix is set at 50 by 50. This 
SUBROUTINE calles to the following subroutines found in SUBS.FOR: 


CSR - Complex square root that preserves negative imaginary 
portion. 

BES1 - Computes the Bessel Function values for a given z. 

HANi - Computes the Hankel Function values for a given z. 


Written 5/88 by R.C. Hurley using MicroSoft Fortran 4.01. 


INTEGER k,m,Mmax,n,Nbig,N1,N2 

REAL l,a,q,d,K2,1T(0:500,0:60),1(¢0:500,1:3),Pi1,R1,R2 
COMPLEX j,zv,Zu,vn,un, JOn, Jin, HOn, Hin, SUM,S$1 
COMPLEX CSR,Amk(Nbig,Nbig), Ja(0:500, 1:4) 

COMPLEX Pmk(0:60,0:500) 

character*1 cr 

COMMON/BLk1/1,j,N1,Pi,d,Nbig 
COMMON/BIk2/K2, | ,Mmax,N2,q,a,T 

cr=char(13) 


This section loads a matrix consisting of Bessel fuctions and 
other values required to load the Amk matrix. The columns are 
filled as follows: 


Ja(n,1) = un*Jo(un*a) in numerator of Sum 

Ja(n,2) = Ji(un*a) in denominator of Sum 

Ja(n,3) = vn*Hi(vn*a) used for all locations 
Ja(n,4) = vn*vn*Ho(vn*a) used in diagonal elements 


WRITE(6,100) Mmax 
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DO 11 n=0,Mmax 
R1=K2-((n*Pi/q)*(n*Pi/q)) 
R2=K2-((n*Pi/l)*(n*Pi/l)) 
un=CSR(R1) 
wn=CSR(R2) 
zu=un*a 
zv=vn*a 
Calt BES1(zu, JOn, Jin) 
Call HAN1(zv,HOn,H1n) 


Ja(n,1) = Un*J0n 

Ja(n,2) = Jin 

Ja(n,3) = Vn*Hin 

Ja(n,4) = Vn*®*Vn*HOn 
c write(7,*) ‘n='’,n,'ja(n,4)=", ja(n,4) 
11 Continue 


WRITE(6,*) ! Completed. / 
WRITEC6,*) ' 
WRITE(6,200) n1,Mmax 


C 
C Loading the Pmk Matrix (size N1 x Mmax) 
6 
DO 22 m0,N1 
WRITE(6,300) m,er 
DO 22 k=0,Mmax 
Pmk(m,k)=(0.0,0.0) 
DO 33 n=0,N2 
S$1=(Ja(n, 1)*T(m,n)*T(k,n))/(CI(n,2)*Ja(n,2)) 
Pmk(m,k)=Pmk(m, k)+S1 
53 CONTINUE 
22 CONT INUE 
WRITEC6,*) ‘ Completed. : 
DO 44 m=0,N1 
DO 44 k=0,N1 


If(m.EQ.k) Then 
Amk (m+ 1, k+1)=Ja(k,3)*Pmk(m,m)-Ja(k,4)*1(m, 1) 


ELSE 
Amk(m+1,k+1)=Ja(k,3)*Pmk(m, k) 
ENDIF 

44 CONTINUE 
100 FORMAT C’ Loading Ja matrix of size’,I13,’ x 4...... ) 
200 = FORMATC! Loading Pmk Matrix of size’,13,’ x‘',13,"..... ) 
300 FORMAT(’ Calculating row’ ,i13,a,\) 

RETURN 

END 
C 
C RAKEKEKEEEEEEEEEEEEEEAKEEEAEEKEEEKKEEEEEEEKEEEEEEEEEEEEEEEKKKKKEKKKKK 
C 

SUBROUTINE BLOAD(FMHZ,M1,Pmk, Ja,B) 
: REEEKEETAKKKKKKKKEKKEKKEKEEKKEKEKKEKEKREEKEKEKKEKEEEREKEEEEEEK 
Cc * & 
c = Variable Definitions z 
Cc * * 
a REEEECEEEEEEEEEEEEEEEETCEEEEETEEEECEEECEEEEECEEEEKEEETEEKETCEKETEKTEEE 
Cc B(m) = Driving Matrix 
Cc N1 = Number of Coefficient Expansions in Region 1 
c I3(m) = Trig. Moment Integrations 
Cc FMHZ = Driving voltage frequency in megahertz 
Cc eps = §8.854E-12 
Cc j = Imaginary One 
C 


Real d,1(0:500,1:3),Pi,eps,w 
Integer k,M,M1,N1,Nbig 
Complex B(Nbig), j,Ja(0:500,1:4),Pmk(0:60,0:500), Sum 
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COMMON/BLK1/1I,j,N1,Pi1,d,Nbig 
eps=8.854E-12 
w=2.0E6*Pi*FMHZ 
DO 22 m-0,N1 
Sum=(0.0,0.0) 
DO 33 k=0,M1 
Sum=Sumt+( (1 (kK, 3)*Jack,3)*Pmk(m,k))/CI Ck, 1)*Jack,4))) 
33 CONTINUE 
B(m+1)=j*w*eps*Sum/d 
22 CONT I NUE 


RETURN 
END 
C 
Cc Pe He Me MH MM HH HH HH HK KKK KKK KEKE KKKKKK 
Cc 
FUNCTION CSR(R) 
REAL R 
COMPLEX CSR 
C Evaluating SQRT(R) with special consideration to neg. reals 
IF(R.GE.0.0) then 
CSR=SQRT(R) 
else 
CSR=(0.,-1.)*SQRT(-R) 
end if 
RETURN 
END 
c 
C Be WH WH KWH HG HH TT HH HE CHK CECE EKKECEKEKKKK 
Cc 
SUBROUTINE BES1(Z,J0,J1) 
c 
c Computing Bessel Functions for n=0,1 with 
C Complex Argument, Z. Direct Power Series Method for 
C CABS(Z) .LE. 6 and Hankel’s Asymptotic Formula for 
Cc CABS(Z) .GT. 6. Written 11/5/87 by M.A. Morgan 
¢ 
INTEGER M,M2 
REAL C(34),D0M, F(34),G0,P(34),Pi,P2 
COMPLEX Z,22,23,24,J0,J1,AM,CL,PO,P1,Q0,Q1,C0,C1,S0,$1 
Pi=3.1415927 
P2=2.0/P! 
IF(CABS(Z).LE.6.0) THEN 
C 
C Utilizing the Direct Power Series Method 
c 
GO= 1.781072 
Z22=0.5*2 
CL=CLOG(G0*Z2) 
Cc 
c Computing F(m) =m! and P(m) = 1 + 1/2 + 1/3 + ....% 1/m 
C 
F(1)=1.0 
P(1)=1.0 
DO 11 M=2,34 
F(M)=M*F(M-1) 
P(M)=P(M-1)+1.0/M 
11 CONTINUE 
C 
c Computing Power Series Coefficients 
c 


DM=-1.0 
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aqn 


qgaagaaAaaAnaAAaA 
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DO 22 M=1,34 
C(M)=DM/(F(M)*F(M)) 
DM=-DM 

CONT I NUE 


Computing JO and J1 


ELSE 


Hankel ! 


ENDIF 


RETURN 


END 


JO=(1.,0.) 
J1=(0.,0.) 

M=0 

M=M+1 

M2=2*M 
AM=C(M)*(Z22**M2) 
J0=J0+AM 
J1=J1-M*AM 

IF ((CABS(AM).GT.1.0E-10).AND.(M.LT.34)) GO TO 33 
J1=J1/Z2 

return 


Asymptotic Formula (Abram. & Stegun p. 364) 


22=2*2Z 

23=2*2Z2 

24=2*23 

P0=1.0- .0703125/22+.1121521/24 
QO0=- .125/2+ .0732422/23 
P1=1.0+.1171875/22- .1441956/24 
Q1=.375/2Z- .10253906/23 
CO=CCOS(Z-.25*PI) 
SO=CSIN(Z-.25*PI) 
C1=CCOS(2-.75*PI) 
S1=CSIN(Z-.75*PI) 
AM=CSQRT(P2/Z) 

JO=AM* (PO*CO0-Q0*S0) 
J1=AM*(P1*C1-Q1*S1) 


RAAEEEAAEEAEAKEEKEEEEEKEAKEEKEKAKEKEKKEEAAAEAEKAKRKEAAKKEEKEEKKKEEKKKEKE 


SUBROUTINE HAN1(Z,HO,H1) 


Computing Hankel Functions for n=0,1 with 

Complex Argument, Z. Direct Power Series Method for 
CABS(Z) .LE. 5 and Hankel’s Asymptotic Formula for 
CABS(Z) .GT. 5. Written 11/6/87 by M.A. Morgan 


INTEGER M,M2 

REAL C(34),DM,F(34),G0,P(34),Pi,P2 

COMPLEX Z,22,23,24,J0,U1,Y0,Y1,AM,CL,P0,P1,Q0,Q1 
COMPLEX E0,E1,X0,X1,HO,H1, j 

P1=3.1415927 

P2=2.0/P1 

j=(0.,1.) 

IF(CABS(Z).LE.5.0) THEN 


Direct Power Series Method 


GO= 1.78072 
22=0.5*2 
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CL=CLOG(G0*Z2) 
C Computing F(m) =m! and P(m) = 1+ 1/2 + 1/3 + ...-+ 1/m 


F(1)=1.0 

P(1)=1.0 

DO 11 M=2,34 
F(M)=M*F(M-1) 
P(M)=P(M-1)+1.0/M 

CONTINUE 


=—_ 


Computing Power Series Coefficients 


qaaono—_ 


DM=-1.0 
DO 22 M=1,34 
C(M)=DM/(F(M)*F(M)) 
DM=-DM 
2 CONT INUE 


Computing JO and J1 


Qa aN 


JO=(1.,0.) 
J1=(0.,0.) 
M=0 
33 M=M+1 
M2=2*M 
AM=C(M)*(2Z2**M2) 
JO=J0+AM 
J1=J1-M*AM 
1F((CABS(AM).GT.1.0E-10).AND.(M.LT.34)) GO TO 33 
J1=J1/22 
C 
C Computing YO and Y1 
Cc 
M=0 
YO=CL* JO 
¥1=Z2*CL*J1-0.5*J0 
44 M=M+1 
M2=2*M 
AM=C(M)*P(M)*(22**M2) 
YO=Y0-AM 
Y1=Y 1+M*AM 
IFC (CABS(AM).GT.1.0E-10).AND.(M.LT.34)) GO TO 44 
YO=P2*Y0 
Y1=P2*Y1/Z2 
HO=J0-j*Y0 
Hi=J1-j*Y1 
RETURN 
ELSE 


C Hankel’ Asymptotic Formula (Abram. & Stegun p. 364 


22=2*2 

23=2*22 

24=2*23 

PO=1.0-.0703125/22+. 1121521/24 
Q0=- . 125/2+.0732422/23 
P1=1.0+.1171875/22- .1441956/24 
Q1=.375/Z- . 10253906/23 
X0=(Z-.25*PI) 

X1=(2-.75*PI) 

EO=CEXP( - j*X0) 

E1=CEXP(- j*X1) 

AM=CSQRT(P2/Z) 
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HO=AM* (PO- j*Q0)*E0 
H1=AM*(P1- j*Q1)*E1 
ENDIF 
RETURN 
END 


C HERRERA EREKREERERKEEKEEEEEEREKREREKREEREREEKEKEEK 


SUBROUTINE FACTOR (A,P,N,NMX) 


Cc PERFORMING LU-DECOMPOSITION WITH PIVOTING ON THE A-MATRIX 
C WRITTEN BY M.A. MORGAN 
C 


INTEGER I,J,JP1,K,P(NMX),PJU,PR,R,RM1,RP1 
REAL DMAX,ELMAG 
COMPLEX AC(NMX,NMX),D(120) 
DO 60 R=1,N 
DO 10 K=1,N 
10 DCK)=A(K,R) 
RM1=R-1 
IFC(RM1.LT.1) GO TO 31 
DO 30 J=1,RM1 
PJ=P(J) 
ACJ,R=DCPU) 
DCPJ)=DCJ) 
JP1=J+1 
DO 20 I=JP1,N 
20 DCI )=DC(I)-ACI, J)*ACU,R) 
30 CONTINUE 
31 CONTINUE 
DMAX=CABS(D(R)) 
PC(R)=R 
RP1=R+1 
IFCRP1.GT.N) GO TO 41 
DO 40 I=RP1,N 
ELMAG=CABS(D(I)) 
IFCELMAG.LT.DMAX) GO TO 40 


DMAX=ELMAG 
40 P(R)=I 
41 CONTINUE 
PR=P(R) 
ACR,R)=D(PR) 
D(PR)=D(R) 


IF(RP1.GT.N) GO TO 51 
DO 50 I=RP1,N 


50 ACI, R)=DC1)/ACR,R) 
51 CONTINUE 
60 CONTINUE 
RETURN 
END 
c 


c RERKKKRAKKKKKKEKEKKEKKKKREKEKKEKREKEKKEKEEKEKKEEEEREREKEEKEKKKEKEKEKKEK 


SUBROUTINE SOLVE (A,B,P,N,NMX) 
c BACK SUBSTITUTION TO INVERT THE LINEAR SYSTEM 
WRITTEN BY M.A. MORGAN 


om 


COMPLEX ACNMX,NMX),8(NMX),Y(120), SUM 
INTEGER I,1P1,U,K,P(NMX),PPI 
DO 20 I=1,N 
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10 


20 


30 
31 


40 


lo] 


MOANA AYAAADANANAANANANANANANANANANANANANAAAAAN 


PPI=P(1) 
Y(1)=B(PP1) 
B(PPI)=B(1) 
IP1=1+1 
IF(IP1.GT.N) GO TO 11 
DO 10 J=IP1,N 
B(J)=B(J)-ACU, 1)*YC(1) 

CONTINUE 

CONTINUE 

CONTINUE 

DO 40 K=1,N 

I=N-K+1 

SUM=(0.,0.) 

IP1=1+1 

IFCIP1.GT.N) GO TO 31 

DO 30 J=IP1,N 
SUM=SUM+A (1, J)*B(J) 

CONTINUE 

CONT INUE 

BCI1)=(Y(1)-SUM)/ACI, 1) 

CONTINUE 

RETURN 

END 


KAAKKAEEKKKKEEEKEKEKKEEKEEKEREKEKEKEEEEKEKKEEEEEEEEEKKKEEEEEEKEKKEEKEEEEEKK 


PROGRAM PLOT 


Program to Plot a Solid Line Plot with Option to Overlay a 
++++ Line Plot for Comparison. Uses PLOT2 Data File Format. 


Original MS-FORTRAN Version of PLOT2 4/24/87 by M.A. Morgan. 
Mods: 5/30/87 XMIN; 11/3/87 RM-FORTRAN 2rd Plot Option. 


INPUT DATA FORMAT 
Unit #3: Solid Line Plot 


TITLE1- 64 Space Header 

N - # Data Points IS (Integer*2) 
XMIN - Real Min X value €£12.3 
XMAX - Real Max X value £E12.3 
FICN) - Input Data Array E12.3 


Unit #4: + Symbol Plot 


TITLE2- 64 Space Header 

N - # Data Points I5 (Integer*2) 
XMIN - Real Min X value €£12.3 
XMAX - Real Max X value £12.3 
F2(N) - Input Data Array €E12.3 


CHARACTER*1 YN,YN1,0UM, YN2, SYMBOL 
CHARACTER*4 LINE 

CHARACTER*7 CROSS 

CHARACTER*16 LTIT,CTIT 
CHARACTER*32 SCALE 

CHARACTER*64 TITLE, FNAME 

REAL X(1025),F1( 1025), F2¢1025) 
INTEGER*2 N,JROW, JCOL, ISYMBL, ITYPE 
SCALE=’is the Vertical Scale Multiplier’ 
LINE='--- ¢ 

CROSS=’ +++ 0! 
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CALL QBEEP 

WRITE(*,*) ‘Press RET to Exit from Screen Plot! 
WRITE(*,*) ‘Printer Hardcopy ? (Y/N): ’ 
READ(*,100) YN 

WRITE(*,*) ‘Enter Data File Name for ---- Plot: ‘ 
READ(*,100) FNAME 

Reading from Unit # 3 

OPEN(3, FILE=FNAME ) 

READ(3,100) TITLE 

READ(3,120) N 

READ(3,130) XMIN 

READ(3,130) XMAX 

WRITE(*,*) ‘---- Plot TITLE, N, XMIN, XMAX: 
WRITE(*,*) TITLE 

WRITEC*,*) N,XMIN, XMAX 

WRITE(*,*) ‘Enter ---- Plot Caption (16 char max): 
READ(*, 100) LTIT 

DX=(XMAX~-XMIN)/(N-171.0) 

FMIN=0.0 

FMAX=0.0 

DO 22 K=1,N 

X(K)=XMIN+(K-1.0)*DX 

READ(3,130) F1(K) 

IFCFICK).LT.FMIN) FMIN=F1(K) 

IFCF1CK).GT.FMAX) FMAX=F1(K) 

F2(K)=0.0 

WRITE(*,*) ‘Comparison Plot (++++) on Same Graph 7??? (Y/N): 
READ(*, 100) YN1 
IFCCYN1.EQ./N’).OR.CYN1.EQ.‘/n’)) GO TO 11 
WRITE(*,*) ‘Enter Data File Name for ++++ Plot: Q 
READ(*, 100) FNAME 

Reading from Unit # 4 

OPEN(4, FILE=FNAME) 

READ(4,100) TITLE 

READ(4,120) N 

READ(4,130) XMIN 

READ(4,130) XMAX 

WRITE(*,*) ‘++++ Plot TITLE, N, XMIN, XMAX:/ 
WRITE(*,*) TITLE 

WRITE(*,*) N,XMIN, XMAX 

WRITE(*,*) ‘Enter ++++ Plot Caption (16 char max): ' 
READ(*,100) CTIT 

DO 27 K=1,N 

READ(4,130) F2(K) 

IFCF2(K).LT.FMIN) FMIN=F2(K) 

IFC F2(K).GT.FMAX) FMAX=F2(K) 

CONTINUE 

CLOSE (4) 

CONTINUE 

CALL QBEEP 

WRITE(*,*) ‘Enter TITLE for Plot (64 Char. Max):/ 
WRITE(*,*) so 

READ (*,100) TITLE 

CLOSE (3) 

IFCFMIN.GT.0.0) FMIN=0.0 

IFC FMAX.LT.0.0) FMAX=0.0 

Computing Scale Factors for Vertical Axis 
ABSMIN=ABS(FMIN) 

ABSMAX=ABS ( FMAX) 

YMAX=AMAX 1 (ABSMIN , ABSMAX ) 
NSCL=INT(LOG10( YMAX) ) 

IF (YMAX.LT.1.0) NSCL=NSCL-1 

YSCL=10.**NSCL 

FMIN=FMIN/YSCL 


oy, 


, 


33 


35 


37 


39 


43 


4&4 


100 
120 
130 
150 


FMAX=FMAX/YSCL 

ABSMIN=ABSMIN/YSCL 
ABSMAX=ABSMAX/YSCL 

DO 33 K=1,N 

F1¢K)=F1ICK)/YSCL 

F2(K)=F2(K)/YSCL 

YMIN=0.0 

IFCFMIN.EQ.0.0) GO TO 37 
YMIN=YMIN+1.0 

IFCABSMIN.GT.YMIN) GO TO 35 
YMIN=YMIN*FMIN/ABSMIN 

CONTINUE 

YMAX=0.0 

IFC FMAX.EQ.0.0) GO TO 41 
YMAX=YMAX+1.0 

IFCABSMAX.GT.YMAX) GO TO 39 
YMAX=YMAX* FMAX/ABSMAX 

CONTINUE 

Calling GRAFMATIC Routines and Plotting F1 Solid Line Graph 
ITYPE=1 

ISYMBL=-2 

CALL QSMODE(6) 

CALL QCMOV(0,21) 

WRITEC*,150) YSCL 

CALL QPTXT(32, SCALE, 3,25, 20) 

CALL QPLOT( 100,530, 12, 147,XMIN, XMAX, YMIN, YMAX,0.,0.,0,1.,1.5) 
CALL QSETUP(O,3,ISYMBL,3) 
XMAJOR=XMAX/5 .0 

CALL QXAXIS(XMIN, XMAX, XMAJOR,1,1,2) 
YMAJOR=1.0 

CALL QYAXISCYMIN, YMAX, YMAJOR,1,1,2) 
CALL QPTXT(64, TITLE,3, 14,24) 

CALL QPTXT(4,LINE,3,14, 22) 

CALL QPTXT(16,LTIT,3, 18,22) 
JROW=(ABSCYMIN)/(CABS( YMAX )+ABS( YMIN) ) )* 135 
JCOL=70-430*XMIN/(XMAX~-XMIN) 

CALL QGTXT(3,/0.0',3,JCOL, JROW,0) 
CALL QTABLCITYPE,N,X,F1) 
IFCCYN1.EQ./N’).OR.CYN1.EQ.’m’)) GO TO 43 
Plotting F2 Graph using + Symbol 
CALL QPTXT(7,CROSS,3,34,22) 

CALL OPTXT(16,CTIT,3,41,22) 
ISYMBL=43 

ITYPE=0 

CALL QSETUP(O,3,ISYMBL,3) 

CALL QTABLCITYPE,N,X,F2) 

CONTINUE 

READ(*,100) DUM 
IFCCYN.EQ./N’).OR.CYN.EQ.’m’)) GO TO 44 
CALL QPSCRN 

CONTINUE 

CALL QSMODE(2) 

WRITEC*,*) ‘End of PLOT Program’ 
FORMAT(A) 

FORMAT(15) 

FORMAT(E12.3) 

FORMAT(13X, 1PE11.3) 

STOP 

END 
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APPENDIX C 


SOFTWARE FLOW CHART FOR MONO WITHOUT GAP FIELD EXPANSION 
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Open 
Fllename 


Write Input 


Varlables to File 


Load Tmn 
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APPENDIX D 


SOURCE CODE FOR MONO WITHOUT GAP FIELD EXPANSION 


PROGRAM MONO 


(Revision 8.1) 


Mod of Mono 7.2, not using the expansion for the Gap Region 
Computing Currents on a Monopole Between Parallel 

Plates using Cylindrical Harmonic Expansions in 2 Regions. 
Links with SUBS4.08J. Created 05/88 LT. R.C. Hurley 


Modified 10/88 M. A. Morgan 


ics oS Sa Varwaltes. ----=------= 


J 
Pi 
KO 
where F 
Cc 
Ke 
N1,N2 
NBIG 
Unit 0 
Unit 6 
Unit 7 


Large: Amk 


Frequency in MHz 

Monopole height in Meters 

Monopole radius in Meters 

Distance between lower ground plane and feed point 
in Meters 

Distance between ground planes in Meters 

height from top of monopole to upper ground plane 
in Meters 

Imaginary One 

Value of pi (3.1415927) 

2*pi*F/c 

frequency of operation in Hz 

Speed of Light 

KO**2 

Number of Coefficients in each of the two regions 
Size of the arrays used in FACTOR and SOLVE 
Standard Input Device 

Screen Output 

Main Output File to Disk (Name selected by User.) 
seaee Variable Definitions ***** 


INTEGER k,m,mpn,n,N1,N2,Nbig,NI ,Nmax,NS,NZ,P(61) 

REAL a,ai,Am,An,c,CI,CM,CP,CR,d,di,DeltaZ,Denom,Dz,Ez2M,eps, FMHZ 
REAL h,hi,1(0:60,1:3),KO,K2,1,li,Lambda,Pi,q,Rin,S 

REAL 1(0:60,0:60),w,Xin,y, yhat,Z,Zmin, Zmax 

COMPLEX Amk(61,61),8(61),C1,C2,1Z(61),j,J0,31,Ja(0:60,1:4) 
COMPLEX Pmk(0:60,0:60),Sum,Y0,Y1,Zin 

CHARACTER*25 Fname 

CHARACTER*1 Bell 

CHARACTER*1 HOW 

CHARACTER*1 INPUT 

CHARACTER*1 Eok 

CHARACTER*12 INname 

COMMON/BLk1/I,j,N1,Pi1,d,Nbig 

COMMON /BLk2/K2,1,Nmax,N2,q,a8,T 


RAEEEKAEKKEREEEEEEEKEEKEREREKEKKKKK 


= Main Program Code . 


RAEAKKERKRERKEKERRERRKREREKEREKREKKKEE 
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REKKAEKKAKKKKAKKKKEKEKEEKEK 


. Input Section sa 
RAEKARAKKKKKKEKKKEKKREKKEKEK 


Bel L=CHAR(7) 
Nbig=61 
c=3.0E+08 
WRITE(6,*) ‘Do you wish to use an input data file [Y/N]?! 
READ(0,100) HOW 
IF (HOW.EQ.’y’ .OR.HOW.EQ.’Y’) THEN 
WRITEC6,*) ‘Enter input data file name with extension :’ 
READ(O,100) INname 
OPEN(Uni t=4, File=INname) 
READ(4,100) INPUT 
IF CCINPUT .EQ. ‘w’) .OR. CINPUT .EQ. ‘’W’) .OR. CINPUT .EQ. 
+ 'm’) .OR. CINPUT .EQ. ’M’)) THEN 
READ(4,*) FMHZ,hi,ai,di, li,N1,N2,NI 
ELSE 
WRITEC*,*) Bell 
WRITE(6,*)’Error detected in first line of data file ’ 
+, INname 
WRITE(6,*) ‘First line must be either M or W for meter 
+ values or wavelength values.’ 
STOP 
ENDIF 
ELSE 
WRITE(6,*) ‘Enter data in either [MJeters or [Wlavelength [M 
+/W)?! 
READ(O,100) INPUT 
IF CINPUT.EQ.’m’ .OR.INPUT.EQ./M’) THEN 
WRITEC6,*) ‘Enter Frequency in MHz: ’ 
READ(O,*) FMHZ 
WRITE(6,*) ‘Enter Monopole Height, h, in meters: 
READ(O,*) hi 
WRITEC6,*) ‘Enter Monopole Radius, in meters: 
READ(O,*) ai 
WRITE(6,*) ‘Enter Feed Gap, d, in meters: 
READ(O,*) di 
WRITEC6,*) ‘Enter Upper Plane Height, | (meters): / 
READ(O,*) Li 
ELSEIF CINPUT.EQ.’w’ .OR.INPUT.EQ./W’) THEN 
WRITE(6,*) ‘Enter Frequency in MHz: ’ 
READ(O,*) FMHZ 
WRITE(6,*) ‘Enter Monopole Height, h, as factor of wav 
+elength : / 
READ(O,*) hi 
WRITE(6,*) ‘Enter Monopole Radius, as factor of wavele 
+ngth : / 
READ(O,*) ai 
WRITE(6,*) ‘Enter Feed Gap, d, as factor of wavelength 


READ(O,*) di 
WRITE(6,*) ‘Enter Upper Plane Height, lt, as factor of 
+wavelength : / 
READ(O,*) ti 
ELSE 
WRITEC*,*) Bell 
WRITEC6,*) ‘You must enter either an M or W. 
+Reenter.. / 
GOTO 2 
ENDIF 
WRITE(6,*) ‘Enter # Region 1 Coeffs. (N1 .LE. 60): d 
READ(O,*) N1 
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WRITE(6,*) ‘Enter # Region 2 Coeffs. (N2 .LE. 60): : 
READ(O,*) N2 
WRITE(6,*) ‘Enter No. Points for I(z) including ends: 
READ(O,*) NI 

ENDIF 


RARAAEAEKAEKEEEEREREEEEKEEEKEKEEKKKKKKKKKKKKKKKKKKK 


* * 
: Calculate the Physical Parameters be! 
* in both Wavelenght and Meters = 
* * 


HACER KERKKKEEREKREREREREREREEREEEKEAKKKKKKKKKKK 


Lambda=c /( FMHZ2* 10**(6)) 
IF CINPUT.EQ.’w’ .OR.INPUT.EQ.’W’) THEN 
h=hi*Lambda 
a=ai*lambda 
d=di*Lambda 
l=li*Lambda 
ELSE 
h=hi 
hi=h/Lambda 
a=al 
aiza/Lambda 
d=di 
di=d/Lambda 
=i 
li=l/Lambda 
ENDIF 


REAAAAEAEAEEEEKEKKKKKEKEKEKEEKEKEERRKREEEKREKEEK 


* * 
: Is Upper Plate Height an se 
ad Even Multiple of Wavelength? se 
* * 


RAEAEAEKKEKAEAAEAKEEKEREEREEAEEEEREKKEEEKEEKKKK 


IF (MOO0((1*2),Lambda).LT. 1£-05) THEN 
WRITEC*,*) Bell 
WRITEC6,*) * ? 
WRITE(6,*) ‘ERROR! Upper Plate Height is almost an even mul 
+tiple of half the wavelength.’ 
WRITE(6,*) ‘Ensure | is NOT within 0.00001 of an integer mul 
+tiple of half the wavelength.’ 
STOP 
ENDIF 


KARA AARKRAKRKKEKKEKKEEEKRERKKKKKKKKKKKKKKKK 


* * 
. Is Upper Plate Height an ss 
= Even Multiple of Antenna Height? * 
* * 


AAARAAEAEKAEKEAEEAEEAEEEEEREEREKEEREEREEKREREEKEE 


IF ((MOO(L,h)).LT.1E-05) THEN 
WRITEC*,*) Bell 
WRITE(6,*) ' ! 
WRITE(6,*) ‘ERROR! Upper Plate Height is almost an even mul 
+tiple of the antenna height.’ 
WRITE(6,*) ‘Ensure lt is NOT within 0.00001 of an integer mul 
+tiple of h (Cin wavelength).’ 
STOP 
ENDIF 
WRITE(6,*) ‘Enter output file name and extension: ' 
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READ(O,100) Fname 

OPEN (Unit=7, Fi le=Fname) 

WRITE(7,210) FMHZ 

WRITE(7,220) h,hi 

WRITE(7,230) a,ai 

WRITE(7,240) d,di 

WRITE(7,250) L,li 

WRITEC(7,*) ‘Region I Coeffs. (N1): ’,N1 
WRITE(7,*) ‘Region II Coeffs. (N2): ’,N2 
WRITE(7,*) ‘No. Points for I(2): NI 
j=(0.0,1.0) 

Pi=3.1415927 

KO=Pi*FMHZ/150.0 

K2=K0*KO 

eps=8 .854E-12 

w=2.0E6*Pi*FMHZ 

yhat=w* eps 

q=l-h 

Nmax=MAX0(N1,N2) 

WRITE(6,*) ‘Loading the Moment Matrices 
WRITE(6,*) ° / 


RAMA RERREAAAAAARAAAAAAAAAAEAAEKEEEKEKAKAKAAEEE 


* 


* Evaluating Moment Integrations 
= Over Region h < z < l. 
* 


* 
* 
* 
* 


Pee eee eee eee eee eee tte eee t eee ee eee eee eee eee eed df | 


T(0,0)=q 
DO 11 k=1,N2 
T(0,k)=0.0 
CONTINUE 
DO 22 m=1,Nmax 
Am=m*Pi/l 
DO 22 n=0,N2 
An=n*Pi/q 


IF(C(ABS(Am-An)).GT.1.0E-3) THEN 


Denom=Am*Am- An*An 

S=SIN(Am*h) 

T(m,n)=-Am*S/Denom 
ELSE 

mon=m+n 


T(m,n)=0.5*q*((-1)**mpn) 


ENDIF 
CONTINUE 


WRITE(6,*) ’ Moment Matrix T Loaded. ’ 


WRITE(6,*) ’ ! 


REMARK AEKAAKKAAKKEKEAKKKAKKREARERARAAKEAKEAEEEKE 


* 


* Loading the Im Matrix for the 
: Largest Number of Coefficients 


® 


* 
* 
* 
* 


Pee eee tee ee ee eee eee eee eee ee eee ee ee ee ee sete ted fo 


00 33 m=0,Nmax 
I1F(m.EQ.0) Then 
I¢(m, 1)=t 
1(m,2)=q 
I1(m,3)=d 
ELSE 
1(m,1)=1/2.0 
I(m,2)=q/2.0 
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33 


40 


Aaaaraan 


59 


aOaaan 


Am=m*Pi/( 
1(m,3)=SINCAm*d)/Am 

ENDIF 
CONTINUE 
WRITE(6,*) / Moment Matrix I Loaded.’ 
WRITE(6,*) ' 
WRITE(6,*) ‘Loading System Matrix .....! 
CALL ALOAD(Ja, Pmk, Amk) 
WRITE(6,*) ‘System Matrix Amk Loaded.’ 
WRITE(6,%) " 
WRITE(6,*) ‘Loading Driving Vector .....! 
DO 40 m0,N1 

B(m+1)=j*yhat*!(m,3)/d 
CONTINUE 
WRITE(6,*) ‘Driving Vector B Loaded.’ 
WRITE(6,*) / # 
WRITE(6,*) ‘Solving the Linear System .. Please Wait .... / 
Ns=N1+1 
CALL FACTOR(Amk,P,Ns,Nbig) 
CALL SOLVE(Amk,B,P,NS,Nbig) 
WRITEC(7,*) ‘Listing the Expansion Coefficients’ 
WRITE(6,*) ¢ Listing the Expansion Coefficients ..... ‘ 
WRITEC7,*) ‘N and BC(N): ° 
DO 44 N=0,N1 

WRITEC7,*) N,BCNt+1) 
CONTINUE 


REAR KKKKKKKKKAKEKKKKKKKKKKEKKEKEEKKKKEKEEEKKEK 


* * 
i Computing Current from the ul 
a Expansion Coefficients id 
® * 


RAEKKKKKKKKKKKKKKKKKEKKKEKEKKKKEEKKKKKEKKKEKEKEKEK 


WRITE(6,*) " * 
WRITE(6,*) ‘Computing current on the monopole..... : 
WRITE(7,*) ‘Listing the Current on the Monopole’ 
DeltaZ=(h-d)/(NI-1) 
WRITE(7,*) ‘Position Distance 1(z) (A)! 
WRITE(7,*) ‘ Number (Meters) Mag Phase’ 
y=d 
k=1 
sum=(0.,0.) 
DO 66 n=0,N1 
Sum=Sum (B(n+1)*Ja(n,3)*cos(n*Pi*y/l)) 
CONTINUE 
1Z2(k)=2.0*P1*a*Sum 
CM=CABS(12(k)) 
CR=REAL(I2Z(k)) 
CI=AIMAG(I2Z(k)) 
CP=180.*ATAN2(CI,CR)/pi 
WRITE(7,200) k,y,CM,CP 
y=y+DeltaZ 
k=k+1 
IF(y.LE.(h+1.E-6)) GO TO 55 
WRITE(6,*) ‘Complete.’ 


RARE EKKEKEKEEKEEKEEEEEEEEEKEEEEKEEKEKEKEE 


* * 
* Computing Input Impedance = 
* * 


RERREKEKEKEKEEEEEEEEEEEEEKEEEEEEEREEEREEKEEEKEKE 


Zin=(1.0,0.0)/12¢1) 


a5 


aOnaoanranannnnanana 


aeOanrandgan 


100 
102 
103 
200 


210 
220 
230 
240 
250 


Rin=REAL(Zin) 

Xin=AIMAG(Zin) 

WRITE(7,*) ‘Input Resistance = ',Rin 
WRITE(7,*) ‘Input Reactance = /,Xin 


REKRAKRREEKREKKEEKRREREEKEKKEKKEEKEEKEKEKKKKEEKREEEK 


* 


* Calculate Electric Field 
* Distribution 
* 


* 
* 
* 
* 


KREKKKEKEKRREEKEKKEKEEEKRERREREERREKEEKKKEEEEEEEREKKEKK 


Checking Ez(a,z) of solution at 200 points from 0 tol 
WRITE(6,*) ‘Do you want to compute the Electric Field (Y/N)? ' 


READ(0,100) Eok 
IF (CEok .EQ. ‘Y’) .OR. (Eok .£Q. ‘y’)) THEN 
NZ=200 
ZMIN=0.0 
ZMAX=L 
D2=Z2MAX/(NZ-1.0) 
OPEN(3, FILE=’EZ.DAT’ ) 
WRITE(3,100) ‘EZ vs. 2’ 
WRITE(3,102) NZ 
WRITE(3,103) ZMIN 
WRITE(3,103) ZMAX 
DO 99 K=1,NZ 
z=(K-1.)*DZ 
c1=(0.,0.) 
DO 88 N=0,N1 


C1=C1+B(N+1)*Ja(N,4)*COS(N*Pi*2/L) 


CONTINUE 

C1=C1/(¢ j*yhat) 

EzM=CABS(C1) 

WRITE(3,103) EzM 
CONTINUE 


WRITE(6,*) ’ Electric Field data stored in file EZ.DAT.’ 


ELSE 
WRITE(6,*) ! Completed. ’ 
ENDIF 


KKKKKEEKEKKKEEKEEKEKKEKEKKKKEKEEKKKEKEEKEEKKEKEKKEKKE 


* 


* 1/O Format Statements 
* 


* 
* 
* 


REKKEKKKKKKEKEKKKKKEKKEEKKEKEKKKEEEEKEKEKKEKEKKKKKK 


FORMAT (A) 
FORMAT(15) 
FORMAT (F 12.3) 


FORMAT (2X,13,5X,F8.3,5X, 1PE12.3,5X,OPF8.2, 1PE12.3,5X, OPF8.2, 


+1PE12.3,5X,OPF8.2) 

FORMAT(’ Driving Frequency, ',9X,F7. 
FORMAT(’ Monopole Height,  ',9X,F7. 
FORMAT(’ Monopole Radius, 1 9X,F 
FORMAT(’ Feed Gap Distance, ’,9X,F 
FORMAT(’ Upper Plane Height,’ ,9X,F 
STOP 

END 


3,’ MHz’) 
7.3,’ meters 
7.3,’ meters 
7.3,’ meters 
7.3,’ meters 
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(’,F5.3,' Lambda)’ ) 
(’,F5.3,' Lambda)’) 
(’,F5.3,’ Lambda)’) 
(’,F5.3,’ Lambda)’) 


qgaaannVannandna 


Subroutine Sorce Code for MONO without Gap Field Expansion 


SUBROUTINE ALOAD(Ja, Pmk, Amk) 


(Revision 8.1) 


RkkkKKKkKkKKkKk&kkkkkk&kkk&kk&kkkhkkkhkhkkk&kkkkkkkkkkk 


* 


* 


* Variable Definitions . 


* 


* 


kkkkkkkk&kkkk&kkkkkkkkkkkkkkkkkkkkkkkkkkkkkhkh 


( = Distance between ground planes 

Radius = Radius of monopole 

q = Distance between top of monopole and upper plane 
d = Antenna feed gap 

K2 = KO**2 (From main program) 

T(m,n) = Moment integration result as function of m and n 
I(m,x) = Moment integration result of order x given m 

j = Imaginary One 

z = Intermediate value for Bessel Subroutine 

un,vwn = Eigenvalues 

JO = Bessel Function of the First Kind of Order 0 
J1 = Bessel Function of the First Kind of Order 1 
HO = Hankel Function of the Second Kind of Order 0 
H1 = Hankel Function of the Second Kind of Order 1 
SUM = Summation over n for the Amk matrix 

$1 = Individual term used in SUM 

m = Column position of matrix Amk 

k = Row position of matrix Amk 

N1 = Number of subdivisons in Region 1 

N2 = Number of subdivisons in Region 2 

Amk = Complex Matrix value that describes the system 


This SUBROUTINE loads the complex A matrix that describes the system 
of size N1 by N1. Maximum size of matrix is set at 50 by 50. This 
SUBROUTINE calles to the following subroutines found in SUBS.FOR: 


CSR - Complex square root that preserves negative imaginary 
portion. 
BES1 - Computes the Bessel Function values for a given z. 


Tmint - Computes the value of the moment integratation T 
Imint - Computes the value of the moment integratation | 
Written 5/88 by R.C. Hurley using MicroSoft Fortran 4.01. 


LARGE: Amk 


INTEGER k,m,n,Nbig,Nmax,N1,N2 

REAL l,a,q,d,K2,1(0:60,0:60),1(¢0:60,1:3),P1,R1,R2 
COMPLEX j,zv,2zu,vn,un, JOn, J1n, HOn, Hin, SUM,S1 
COMPLEX CSR,Amk(Nbig,Nbig),Ja(0:60,1:4) 


COMPLEX Pmk(0:60,0:60) 


CHARACTER*1 cr 


COMMON/BLK1/1,j,N1,Pi,d,Nbig 
COMMON/BLk2/K2,1,Nmax,N2,q,a,T 


cr=CHAR( 13) 


This section loads a matrix consisting of Bessel fuctions and 
other values required to load the Amk matrix. The columns are 


filled as follows: 


Ja(n, 1) 
Ja(n,2) 
Ja(n,3) 
Ja(n,4) 


un* Jo(un*a) in numerator of Sum 

J1(un*a) in denominator of Sum 
vn*Hi(vn*a) used for all locations 
vn*vn*Ho(vn*a) used in diagonal elements 
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WRITE(6,100) Nmax 

DO 11 n=0,Nmax 
R1=K2-((n*Pi/q)* (n*Pi/q)) 
R2=K2-((n*Pi/1)*(n*Pi/L)) 
un=CSR(R1) 
vwn=CSR(R2) 
ZU=un* 8 
ZV=vn*a 
Call BES1(zu, JOn, JI1n) 
Call HAN1(zv,HOn,H1n) 
Ja(n,1) = Un*JO0n 


Ja(n,2) = Jin 

Ja(n,3) = Vn*Hin 

Ja(n,4) = Vn*Vn*HOn 
11 Cont inue 


WRITE(6,*) ' Completed. ’ 
WRITE(6,*) * # 
WRITE(6,200) N1,Nmax 


G 
C Loading the Pmk Matrix (size N1 x Nmax) 
c 
DO 22 m=0,N1 
WRITE(6,300) m,cr 
DO 22 k=0,Nmax 
Pmk(m,k)=(0.0;0.0) 
DO 33 n=0,N2 
$1=(Ja(n, 1)*T¢(m,n)*T¢(k,n))/CI(n, 2)*Ja(n,2)) 
Pmk (m, k )=Pmk(m, k)+S1 
33 CONTINUE 
22 CONTINUE 
WRITE(6,*) ' Completed. : 
DO 44 m=0,N1 
DO 44 k=0,N1 


If(m.EQ.k) Then 
Amk (m+1,k+1)=Ja¢(k,3)*Pmk(m,m)-Jack,4)*1¢m, 1) 


ELSE 
Amk (m1, k+1)=Jack,3)*Pmk(m,k) 
ENDIF 
44 CONTINUE 
100 FORMAT(‘ Loading Ja matrix of size’,13,’ x 4...... r) 
200 FORMAT(‘ Loading Pmk Matrix of size’,13,/’ x',13,’...+-') 
300 FORMAT(’ Calculating row’,13,8,\) 
RETURN 
END 
C 
C KREKKRKREKEKEKKKEKEKEKKKKKKRKKEKRKKKEEKKKKRKKKKRKKKKRKKEKKKKKKKKRKKKKKKKKKKKKKK 
C 
FUNCTION CSRC(R) 
REAL R 
COMPLEX CSR 
c Evaluating SQRT(R) with special consideration to neg. reals 


IF(R.GE.0.0) THEN 
CSR=SQRT(R) 
ELSE 
CSR=(0.,-1.)*SQRT(-R) 
END IF 
RETURN 
END 
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REEEEKKKEKKEKKKKEKEEEKEKEKKKKEKEKKKEKKEEEEEEKEKKKEEKKEEKEKEKEEEREKREKEKE 


aon 


SUBROUTINE BES1(Z, JO, J1) 


Computing Bessel Functions for n=0,1 with 

Complex Argument, Z. Direct Power Series Method for 
CABS(Z) .LE. 6 and Hankel’s Asymptotic Formula for 
CABS(Z) .GT. 6. Written 11/5/87 by M.A. Morgan 


qganaay»n 


INTEGER M,M2 

REAL C(34),DM,F(34),G0,P(34),Pi,P2 

COMPLEX 2,22,23,24,J0,J1,AM,CL,PO,P1,Q0,Q1,C0,C1,S0,$1 
Pi=3.1415927 

P2=2.0/PI 

IF(CABS(Z).LE.6.0) THEN 


wo) 


Utilizing the Direct Power Series Method 


GO= 1.781072 
22=0.5*2 
CL=CLOG(G0*22) 


ie] 


Computing F(m) =m! and P(m) = 1+ 1/2 + 1/3 + ....+ 1/m 


F(1)=1.0 

P(1)=1.0 

DO 11 M=2,34 
F(M)=M*FC(M-1) 
P(M)=P(M-1)+1.0/M 

CONTINUE 


he 


Computing Power Series Coefficients 


qaaaj— 


DM=-1.0 
DO 22 M=1,34 
C(M)=DM/(F(M)*FCM)) 
DM=-DM 
2 CONTINUE 


Computing JO and J1 


aa anv 


JO=(1.,0.) 
J1=(0.,0.) 
M=0 

33 M=M+1 
M2=2*M 
AM=C(M)*(22**M2) 
JO=J0+AM 
J1=J1-M*AM 
1F((CABS(AM).GT.1.0E-10).AND.(M.LT.34)) GO TO 33 
J1=J1/22 
return 

ELSE 


C Hankel’ Asymptotic Formula (Abram. & Stegun p. 364) 


Z22=2*2 

23=2*22 

24=2*23 

PO0=1.0- .0703125/22+ .1121521/24 
Q0=- .125/2+ .0732422/23 


Slee, 


P1=1.0+.1171875/22- .1441956/24 
Q1=.375/Z- .10253906/23 
CO=CCOS(2-.25*PI) 
SO=CSIN(Z-.25*P1) 
C1=CCOS(Z-./75*PI) 
S1=CSIN(2-.75*PI1) 
AM=CSOQRT(P2/Z) 
JO=AM* (PO*CO-Q0*S0) 
J1=AM* (P1*C1-01*S1) 

ENDIF 

RETURN 

END 


Cc * RHREARKREAEEEKRKKEKAEEAEKEKKEEAEEREKKKEEEEEKKERKEKEEKKKEKEEKEKKEKKKEEK 


c 
SUBROUTINE HAN1(Z,HO,H1) 


Computing Hankel Functions for n=0,1 with 

Complex Argument, Z. Direct Power Series Method for 
CABS(Z) .LE. 5 and Hankel’s Asymptotic Formula for 
CABS(Z) .GT. 5. Written 11/6/87 by M.A. Morgan 


angrAaAaANnA 


INTEGER M,M2 

REAL C(34),DM,F(34),GO0,P(34),Pi,P2 

COMPLEX Z,22,23,24,J0,J1,Y0,Y1,AM,CL,P0,P1,Q0,Q1 
COMPLEX E0,E1,X0,X1,HO,H1, j 

P1=3. 1415927 

P2=2.0/P1 

j=(0.,1.) 

IF(CABS(Z).LE.5.0) THEN 


o>] 


Direct Power Series Method 


GO= 1.78072 
22=0.5*Z 
CL=CLOG(G0*Z2) 


] 


Computing F(m) = m! and P(m) = 1 + 1/2 + 1/3 + ....+ 1/m 


F(1)=1.0 
P(1)=1.0 
DO 11 M=2,34 
F(M)=M*F(M-1) 
P(M)=P(M-1)+1.0/M 
1 CONTINUE 


Computing Power Series Coefficients 


aoan— 


DM=-1.0 
DO 22 M=1,34 
C(M)=DM/(F(M)*F(M)) 
DM=-DM 
2 CONT INUE 


Computing JO and J1 


aanrIn 


JO=(1.,0.) 
J1=(0.,0.) 
M=0 

33 M=M+1 
M2=2*M 
AM=C(M)*(22**M2) 


100 


C 
C 
C 


4s 


Cc 


10 


JO=JO+AM 

J1=J1-M*AM 
IF(CCABS(AM).GT.1.0E-10).AND.(M.LT.34)) GO TO 33 
J1=J1/Z2 


Computing YO and Y1 


M=0 
YO=CL* JO 
Y1=Z2*CL*J1-0.5*J0 
M=M+1 
M2=2*M 
AM=C(M)*P(M)*(Z2**M2) 
YO=Y0-AM 
Y1=Y1+M*AM 
IF ((CABSCAM).GT.1.0E-10).AND.(M.LT.34)) GO TO 44 
YO=P2*Y0 
Y1=P2*Y1/2Z2 
HO=J0- j*Y0 
H1i=J1-j*Y1 
RETURN 

ELSE 


Hankel’ Asymptotic Formula (Abram. & Stegun p. 364 


22=2*Z 
23=2*Z2 
24=2*23 
PO=1.0- .0703125/Z2+.1121521/24 
Q0=- .125/2+ .0732422/23 
P1=1.0+.1171875/Z2- .1441956/24 
Q1=.375/2- .10253906/23 
X0=(Z- .25*PI) 
X1=(Z-.75*PI) 
EQ=CEXP(- j*X0) 
E1=CEXP(- j*X1) 
AM=CSQRT(P2/Z) 
HO=AM* (PO- §*Q0)*E0 
H1=AM* (P1- §*Q1)*E1 

ENDIF 

RETURN 

END 


KAAKEKEKKERKKEKKKEKKEKEKAKKKKAKEEKKKEKEEKEKEKEKEKEEEKEKEKREREKKKKKEEKK 


SUBROUTINE FACTOR (A,P,N,NMX) 
PERFORMING LU-DECOMPOSITION WITH PIVOTING ON THE A-MATRIX 


INTEGER I,J,JP1,K,PCNMX),PJ,PR,R,RM1,RP1 
REAL DMAX,ELMAG 
COMPLEX A(NMX,NMX),D(120) 
DO 60 R=1,N 
DO 10 K=1,N 
O(K)=A(K,R) 
RM1=R-1 
IF(RM1.LT.1) GO TO 31 
DO 30 J=1,RM1 
PJ=P(J) 
ACJ,R)=D(Pd) 
D(PJ)=D(J) 
JP1=J+1 
DO 20 I=uP1,N 
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20 


31 


40 
41 


50 


60 


C 
C 
Cc 


10 


20 


30 
31 


40 


D(T)=D(1)-ACI, J) *ACJ,R) 
CONTINUE 
CONTINUE 
DMAX=CABS(D(R)) 
P(R)=R 
RP1=R+1 
IFCRP1.GT.N) GO TO 41 
DO 40 I=RP1,N 
ELMAG=CABS(D(1I)) 
IF(ELMAG.LT.DMAX) GO TO 40 
DMAX=ELMAG 
P(R)=1 
CONTINUE 
PR=P(R) 
ACR,R)=D(PR) 
D(PR)=D(R) 
IFCRP1.GT.N) GO TO 51 
DO 50 I=RP1,N 
ACT,R)=D(1)/ACR,R) 
CONTINUE 
CONTINUE 
RETURN 
END 


KARE KKKKKEKKKKKK KK KKK KKK KKK KK KKK KKEKEKEKCKKKCKCKKKCKKKKKKEKKKEKKKKEK 


SUBROUTINE SOLVE (A,B,P,N,NMX) 
BACK SUBSTITUTION TO INVERT THE LINEAR SYSTEM 
COMPLEX A(NMX,NMX),BCNMX),Y(120),SUM 
INTEGER I,1P1,J,K,PC(NMX),PPI 
DO 20 I=1,N 
PPI=P(I) 
YCI)=BCPPI) 
BC PPI )=BC(1) 
IP1l=I+1 
IFCIP1.GT.N) GO TO 11 
DO 10 J=IP1,N 
BC J)=BC(J)-ACJ, 1) *YC1) 
CONTINUE 
CONTINUE 
CONT INUE 
DO 40 K=1,N 
I=N-K+1 
SUM=(0.,0.) 
IP1=I+1 
IFCIP1.GT.N) GO TO 31 
DO 30 J=IP1,N 
SUM=SUM+A(I, J)*BCJ) 
CONTINUE 
CONTINUE 
BCI )=(YC(1)-SUM)/ACI,1) 
CONTINUE 
RETURN 
END 
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APPENDIX E 


SAMPLE INPUT/OUTPUT DATA FILES USED WITH MONO 


This Appendix lists the input and output data file 
structures used within MONO. There are two different types 
of input data files, one that is entered at the keyboard at 
run time and one which allows unattended processing of a 
number of runs without user intervention. In each case, the 
data fields may be either arranged in column format with one 
entry per line, or all on one line separated by commas. 
Additionally, a variable type is set by the program which 
allows a real variable to read both an input of "1" or "1." 


as the same. 


A. INPUT DATA FILE STRUCTURE 


This data file may be adjusted to allow the user the 
option of selecting either the output file name or, for the 
case of MONO8, the calculation of the electric field along the 
antenna at run time. Below is a commented data set that could 
be used. Note the comments are NOT to be entered in the data 


set as errors would occur. 
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Data File 
Entry 


W 


299m 


~24 


~O1 
-02 
Zoe 
30 
30 


60 


30 


d:\file.out 


eve or uw n" 


Comment 


This identifies the basis for all physical 
measurements of the antenna. Using a "w" 
or "W" signifies all measurements are in 
terms of wavelength, while a "m" or "M" 
indicates measurements are in terms of 
meters. 


Driving source frequency in MHz 


Antenna height, h. Can be entered as 
WO..24 


Antenna radius, a. 

Gap distance, q. 

Ground plate spacing, /. 

Number of modes for Region I, N,. 
Number of modes for Region II, N,. 


Number of modes for Gap expansion, M,. 
Not used for MONO8. 


Number of points to calculate the current 
along the antenna. 


DOS path and name to assign to output 
file. Maximum length is 25 characters. 
This entry may also be entered at the key- 
board at run time. 


Indicates whether to calculate the elec- 
tric field along the antenna. Test is NOT 
case sensitive. If yes, data is stored as 
'EZ.DAT' in the current directory in for- 
mat for PLOT.FOR found in Appendix B. 

This may be entered via the keyboard at 
run time. 


The same data file can be arranged as follows: 


w,299.8,.24,.01.02,2.17,30,30,60,30,d:\file.out,n 
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The data file is easily modified to allow the user to 
group a series of runs together to be batch processed later 
using the standard DOS redirection command. For example, to 
use input data file named 'RUN1' and to have the normal screen 
output redirected to a running log file named 'OUTPUT.LOG', 
the correct syntax would be: 

MONO <[path]RUN1 >>[path]OUTPUT. LOG 

where path 1s the DOS path where the individual files are 
stored. This parameter is optional and defaults to the 
current directory. A number of runs could then be completed 
without user input by combining a series of these commands in 
one batch file. The file structure is the same as that listed 
above except that an 'n' is inserted at the beginning of the 
file to indicate that a separate input file is not to be used. 
Therefore, the same data file used above would be modified as 
seen below: 


nee OOM nOn O22 nl 7cO80,60,30,d:\file.out,n 


B. SAMPLE OUTPUT FILE 


The output file generated by MONO7 for the inputs used 
above is listed below: (* indicates not included in output 


from MONO8) 


ckOD 


Driving Frequency, 299.800 MHz 


Monopole Height, 
Monopole Radius, 


.240 meters 
-010 meters 


Ci Et TE tt all) 


Feed Gap Distance, .010 meters 
Upper Plane Height, 2.171 meters 
Region I Coeffs. (N1): 30 
Region II Coeffs. (N2): 30 
No. Gap Coeffs. (M1): 60 
No. Points for I(z): 30 


Listing the Expansion Coefficients 


N and B(N): 


WOON AO FWNF O 


30 


(-9.259474E-05,-3.712934E-04) 
(-1.716743E-04,-7.285348E-04) 
(-1.250944E-04,-6.843712E-04) 
(-8.932992E-06,-5.970440E-04) 
(6.372479E-04, -3.034107E-04) 
(-7.747178E-04,-5.976821E-04) 
(-4.176759E-04,-5.106181E-04) 
(-2.877242E-04, -4.189367E-04) 
(-2.029108E-04, -3.268855E-04) 
(-1.371401E-04,-2.384746E-04) 
(-8.359972E-05, -1.573042E-04) 
(-4.017067E-05,-8.621241E-05) 
(-6.131331E-06,-2.729695E-05) 

(1.893293E-05,1.824955E-05) 

(3.552895E-05,5.013411E-05) 

(4.440288E-05,6.891159E-05) 

(4.656793E-05,7.587663E-05) 

(4.326494E-05,7.290936E-05) 

(3.588156E-05, 6. 228465E-05) 

(2.588031E-05, 4.649168E-05) 

(1.465865E-05, 2.799412E-05) 

(3.497886E-06,9.087322E-06) 
(-6.533222E-06,-8.273393E-06) 
(-1.463069E-05, -2.258123E-05) 
(-2.028524E-05, -3.284282E-05) 
(-2.328759E-05,-3.859688E-05) 
(-2.370753E-05, -3.988246E-05) 
(-2.185403E-05,-3.717663E-05) 
(-1.821469E-05,-3.127903E-05) 
(-1.339421E-05,-2.322528E-05) 
(-8.037723E-06,-1.412969E-05) 


Listing the Current on the Monopole 


Position 
Number 


OOF WN 


-240 Lambda) 
.010 Lambda) 
-010 Lambda) 


2.170 Lambda) 


Distance I(z) (A) 
(Meters) Mag Phase 
-010 2.316E-02 -15.04 
S018 2. 309E-02 -15.98 
-026 2.299E-02 -17.24 
~034 2.284E-02 -18.58 
-042 2.262E-02 -19.80 
-050 2.234E-02 -20.78 
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7 ~-058 2.200E-02 Sol. 49 


8 066 2.161E-02 -21.97 
9 .073 2.121E-02 -22.34 
10 -081 2.082E-02 22 2 
11 -089 2.044E-02 S86 ALE, 
12 -097 2.007E-02 235.10 
13 -105 1.970E-02 -24.44 
14 i lie gs 1.929E=-02 = 2/5 ol 
15 ee | 1.881E-02 =~29-70 
16 7129 1.822E-02 = 216 506 
17, ~137 Teo O2E=O2 -26.50 
18 ~145 1.669E-02 -26.74 
19 e153 1.574E-02 iow oD 
20 -161 1.467E-02 =U ANS, 
21 - 169 1.349E-02 -27.49 
22 ~177 1.222E-02 -27.86 
23 -185 1.088E-02 2 Oma 
24 ~193 9.485E-03 -28.64 
25 -200 8.069E-03 -28.91 
26 -208 6.671E-03 -29.03 
27 Z2LG 5.331E-03 =28.99 
28 ~224 4.091E-03 -28.84 
29 e232 2.987E-03 =28.72 
30 -240 2.042E-03 -28.84 
Input Resistance = 41.703650 
Input Reactance = 11.206850 
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APPENDIX F 


SAMPLE INPUT/OUTPUT DATA FILES USED WITH NEC 


This Appendix lists an input data file used NEC and the 
corresponding output file. The run selected was for an 


antenna with the following physical characteristics using one 


reflection: 
Driving Frequency 299.8 MHz 
Antenna Height 0.48 wavelengths 
Antenna Radius 0.01 wavelengths 
Gap Distance 0.02 wavelenghts 
Ground Plate Spacing 1.4 wavelengths 

CE Run 11: Monopole Antenna of Height .48,Radius .01,.02,1.40 


GW 1, 24, 0, 0,0, 0, Ceaser 

GW 2, 48, 0, OF) 2.32, Ope ee 2 oyeeeol 
GE 1 

GN 1 

EK 


NO 


0 
EX 0 24 
EX 0 
XQ 
EN 


cans 
L] 

Rpt 

Sree: 
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Sample output file from the above data is shown below: 


RARER KKRERKKKKKRKREREKKEKKKKKE 


NUMERICAL ELECTROMAGNETICS CODE 


KARR KAKA KKKKKRAKRKKKAKRKAKKKKKK 


- - - - COMMENTS - - - - 


Run 11: Monopole Antenna of Height .48,Radius .01,.02,1.40 


== > STRUCTURE SPECIFICATION - ->- 


COORDINATES MUST BE INPUT IN 
METERS OR BE SCALED TO METERS 
BEFORE STRUCTURE INPUT IS ENDED 


WIRE NO. OF FIRST 
LAST TAG 
NO. X1 rl Z1 X2 2 Z2 RADIUS SEG. SEG. 
SEG. NO. 
ak .00000 .00000 .00000 .00000 .00000 - 48000 .01000 24 1 
24 Bl 
2 .00000 .00000 2.32000 .00000 .00000 3.28000 .01000 48 25 
72 Z 


GROUND PLANE SPECIFIED. 


WHERE WIRE ENDS TOUCH GROUND, CURRENT WILL BE INTERPOLATED TO IMAGE IN GROUND PLANE. 


TOTAL SEGMENTS USED= 72 NO. SEG. IN A SYMMETRIC CELL= 72 SYMMETRY FLAG= 0 


- MULTIPLE WIRE JUNCTIONS - 
JUNCTION SEGMENTS (- FOR END 1, + FOR END 2) 


NONE 
Suse = SEGMENTATION DATA - - - = 
COORDINATES IN METERS 
I+ AND I- INDICATE THE SEGMENTSBEFORE AND AFTER I 
SEG. COORDINATES OF SEG. CENTER SEG. ORIENTATION ANGLES WIRE CONNECTION DATA TAG 
NO. X ¥ Z LENGTH ALPHA BETA RADIUS i I Is NO. 
ul .00000 .00000 .01000 .02000 90.00000 .00000 .01000 1 1 Z uN 
2 .00000 .00000 .03000 .02000 90.00000 .00000 .01000 1 Zz 3 nH 
3 .00000 .00000 .05000 .02000 90.00000 .00000 .01000 2 3 4 i 
4 .00000 .00000 .07000 .02000 90.00000 .00000 .01000 3 4 5 u 
> .00000 .00000 .09000 .02000 90.00000 00000 .01000 & 3 6 1 
6 .00000 .00000 . 11000 .02000 90.00000 00000 .01000 be 6 7 1 
7 .00000 .00000 . 13000 .02000 90.00000 .00000 .01000 6 7 8 1 
8 .00000 .00000 . 15000 .02000 90.00000 .00000 .01000 7 8 9 i 
se) .00000 .00000 - 17000 .02000 90.00000 .00000 .01000 8 SS) 10 1 
10 .00000 .90000 . 19000 .02000 90.00000 .00000 .01000 9g 10 11 1 
11 .00000 .00000 . 21000 .02000 90.00000 .00000 .01000 10 11 Vz H 
12 .00000 .00000 . 23000 .02000 90.00000 .00000 .91000 11 12 13 1 
13 .00000 .00000 . 25000 .02000 90.00000 .00000 .01000 12 13 14 al 
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14 .00000 .00000 
NS) .00000 .00000 
16 .00000 .00000 
17 .00000 .00000 
18 .00000 .00000 
19 . 00000 .00000 
20 .00000 .00000 
21 .00000 .00000 
22 .00000 .00000 
23 .00000 .00000 
24 .00000 .00000 
Zo . 00000 . 00000 
26 .00000 .00000 
27 .00000 .00000 
28 .00000 .00000 
29 .00000 .00000 
30 .00000 .00000 
31 .00000 .00000 
32 .00000 .00000 
33 .00000 .00000 
34 .00000 .00000 
35 .00000 .00000 
36 .00000 .00000 
37 .00000 .00000 
38 .00000 .00000 
39 .00000 .00000 
40 .00000 .00000 
41 .00000 .00000 
42 .00000 . 00000 
43 .00000 .00000 
44 .00000 .00000 
45 .00000 .00000 
46 .00000 .00000 
47 .00000 .00000 
48 .00000 .00000 
49 .00000 .00000 
50 .00000 .00000 
ot .00000 .00000 
af .00000 .00000 
53 .00000 . 00000 
54 .00000 .00000 
ao .00000 .00000 
56 .00000 .00000 
57 .00000 .00000 
58 .00000 .00000 
59 .00000 .00000 
60 .00000 .00000 
61 .00000 .00000 
62 .00000 .00000 
63 .00000 .00000 
64 .00000 .00000 
65 .00000 .00000 
66 .00000 .00000 
67 .00000 .00000 
68 .00000 .00000 
69 .00000 .00000 
70 .00000 .00000 
71 .00000 .00000 
72 .00000 .00000 


wex** DATA GARD NO. 1 
0.00000E+00 0.00000E+00 
**xx*** DATA CARD NO. 2 
0.00000E+00 0.00000E+00 
*x*x** DATA CARD NO. 3 
0.00000E+00 0.00000E+00 
***x** DATA CARD NO. 4 
0.00000E+00 0.00000E+00 
*x*e* DATA CARD NO. 5 
0.00000E+00 0.00000E+00 


RRR RR 


OwWWWVUYWWWWWWWWHWNHNNNNNNNNNNNNNNNNNNNNNHNNNNYNNN NNN N PY 


. 27000 
. 29000 
. 31000 
. 33000 
. 35000 
. 37000 
. 39000 
41000 
- 43000 
45000 
-47000 
. 33000 
. 35000 
.37000 
. 39000 
.41000 
43000 
-45000 
. 47000 
. 49000 
. 91000 
. 93000 
- 55000 
. 57000 
. 99000 
.61000 
63000 
.65000 
.67000 
.69000 
. 71000 
. 73000 
. 75000 
.77000 
. 79000 
.81000 
.83000 
85000 
. 87000 
.89000 
.91000 
.93000 
. 95000 
.97000 
. 39000 
.01000 
.03000 
.05000 
.07000 
.09000 
. 11000 
. 13000 
. 15000 
.17000 
. 19000 
.21000 
. 23000 
.25000 
. 27000 


.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 
.02000 


24 
22 


90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 


90. 


90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 
90 


0 0.00000E+00 
0 0.00000E+00 
0 1.00000E+00 
0 -1.00000E+00 


0 -1.00000E+00 


.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
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.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.90000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 
.00000 


.01000 13 14 
.01000 14 15 
.01000 iS 16 
.01000 16 17 
.01000 17 18 
.01000 18 19 
.01000 19 20 
.01000 20 21 
.01000 21 22 
.01000 ZZ 23 
.01000 23 24 
.01000 0 Zo 
.01000 25 26 
.01000 26 27 
.01000 27 28 
.01000 28 29 
.01000 29 30 
.01000 30 31 
.01000 31 32 
.01000 32 33 
.01000 33 34 
.01000 34 35 
.01000 35 36 
.01000 36 37 
.01000 37 38 
.01000 38 39 
.01000 39 40 
.01000 40 41 
.01000 41 42 
.01000 42 43 
.01000 43 44 
.01000 44 45 
.01000 45 46 
.01000 46 47 
.01000 47 48 
.01000 48 49 
.01000 49 50 
.01000 50 spl 
.01000 a1 32 
.01000 a2 53 
.01000 53 54 
.01000 54 oa 
.01000 5)5) 56 
.01000 56 57 
.01000 a7 58 
.01000 58 39 
.01000 39 60 
.01000 60 61 
.01000 61 62 
.01000 62 63 
.01000 63 64 
.01000 64 65 
.01000 65 66 
.01000 66 67 
.01000 67 68 
.01000 68 69 
.01000 69 70 
.01000 70 71 
.01000 71 72 


0.00000E+00 0.00000E+00 


0.00000E+00 0.00000E+00 


0.00000E+00 0.00000E+00 


0.Q0000E+00 0.00000E+00 


0.00000E+00 0.00000E+00 


= 
tA 
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNP PRP RRR PrP re 


.00000E+00 
. O0000E+00 
.00000E+00 
. 00000E+00 


. 00000E+00 


waxee DATA CARD NO. 6 XQ O 0 0 0 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 
0.Q00000E+00 0.00000E+00 


FREQUENCY= 2.9980E+02 MHZ 
WAVELENGTH= 1.0000E+00 METERS 


APPROXIMATE INTEGRATION EMPLOYED FOR SEGMENTS MORE THAN 1.000 WAVELENGTHS APART 
THE EXTENDED THIN WIRE KERNEL WILL BE USED 


- - - STRUCTURE IMPEDANCE LOADING - - - 


THIS STRUCTURE IS NOT LOADED 


- - - ANTENNA ENVIRONMENT - - - 


PERFECT GROUND 


- - - MATRIX TIMING - - - 


FILL= 1.500 MIN., FACTOR= .364 MIN. 


- - - ANTENNA INPUT PARAMETERS - - - 


TAG SEG. VOLTAGE (VOLTS) CURRENT (AMPS) IMPEDANCE (OHMS) ADMITTANCE 
(MHOS ) POWER 
NO. NO. REAL IMAG. REAL IMAG. REAL IMAG. REAL 
IMAG. (WATTS) 
1 1 1.00000E+00 0.00000E+00 2.43054E-03 2.78753E-03 1.77699E+02-2.03799E+02 2.43054E-03 
2.787 53E-O37.215Z75-03 
2 48-1.Q0000E+00 0.00000E+00-2.42866E-03-2.79143E-03 1.77398E+02-2.03896E+02 2.42866E-03 
2.79143E-03 1.21433E-03 
2 49-1.00000E+00 0.00000E+00-2.42265E-03-2.79230E-03 1.77274E+02-2.04322E+02 2.42265E-03 


2.79230E-03 1.21132E-03 , 


- - - CURRENTS AND LOCATION - - - 


DISTANCES IN WAVELENGTHS 


SEG. TAG COORD. OF SEG. CENTER SEG. - - - CURRENT (AMPS) - - - 
NO. NO. Xx y Z LENGTH REAL IMAG. MAG. PHASE 
I a .0000 .0000 .0100 .02000 2.4305E-03 2.7875E-03 3.6984E-03 48.914 
2 a .0000 -0000 .0300 .02000 2.4147E-03 1.0516E-03 2.6338E-03 23.534 
3 a .0000 .0000 .0500 -02000 2-3833E-03 -8.5797E-05 2.3848E-03 -2.062 
& 1 .0000 .0000 .0700 .02000 2.3366E-03 -9.8162E-04 2.5344E-03 ~-22.787 
5 1 .0000 .0000 .0900 .02000 2.2754E-03 -1.7685E-03 2.8819E-03 -37.856 
6 1 .0000 .0000 .1100 .02000 2.2004E-03 -2.4668E-03 3.3055E-03 -48.267 
7 1 .0000 .0000 . 1300 .02000 2.1125E-03 -3.0847E-03 3.7387E-03 -55.595 
8 a .0000 .0000 . 1500 .02000 2.0130E-03 -3.6240E-03 4.1456E-03 -60.950 
9 1 .0000 .0000 .1700 .02000 1.9031E-03 -4.0841E-03 4.5058E-03 -65.016 
10 1 .0000 .0000 . 1900 .02000 1.7842E-03 -4.4630E-03 4.8065E-03 ~-68.209 
11 1 .0000 .0000 .2100 .02000 1.6579E-03 -4.7589E-03 5.0394E-03 -70.792 
12 1 .0000 .0000 .2300 .02000 1.5258E-03 -4.9697E-03 5.1987E-03 -72.932 
13 1 .0000 .0000 .2500 .02000 1.3895E-03 -5.0945E-03 5.2806E-03 -74.744 
14 1 .0000 .0000 .2700 .02000 1.2508E-03 -5.1325E-03 5.2827E-03 -76.304 
is 1 .0000 .0000 . 2900 .02000 1.1113E-03 -5.0841E-03 5§.2041E-03 -77.670 
16 1 .0000 .0000 .3100 .02000 9.7278E-04 -4.9502E-03 5.0449E-03 -78.882 


lil 


17 1 .0000 .0000 . 3300 .02000 8.3692E-04 -4.7327E-03 4.8061E-03 -79.972 
18 J .0000 .0000 . 3500 .02000 7.0529E-04 -4.4340E-03 4.4897E-03 -80.962 
19 at .0000 .0000 .3700 .02000 5.7941E-04 -4,0569E-03 4.0981E-03 -81.872 
20 1 .0000 .0000 .3900 .02000 4.6064E-04 -3.6044E-03 3.6337E-03 -82.717 
21 i .0000 .0000 . 4100 .02000 3.5018E-04 -3.0785E-03 3.0983E-03 <-83.510 
22 i .0000 .0000 . 4300 .02000 2.4894E-04 -2.4785E-03 2.4909E-03 -84.264 
23 1 .0000 .0000 4500 .02000 1.5733E-04 -1.7963E-03 1.8032E-03 -84.995 
24 2 .0000 .0000 .4700 .02000 7.2003E-05 -9.5821E-04 9.6091E-04 -85.703 
25 2 .0000 .0000 2.3300 .02000 -8.2718E-05 9.5969E-04 9.6325E-04 94.926 
26 2 .0000 .0000 2.3500 .02000 -1.7730E-04 1.7990E-03 1.8077E-03 95.629 
27 2 .0000 .0000 2.3700 .02000 -2.7635E-04 2.4819E-03 2.4972E-03 96.353 
28 2 .0000 .0000 2.3900 .02000 -3.8402E-04 3.0826E-03 3.1064E-03 97.101 
29 2 .0000 .0000 2.4100 .02000 ~-5.0002E-04 3.6090E-03 3.6434E-03 97.888 
30 2 .0000 .0000 2.4300 .02000 -6.2347E-04 4.0618E-03 4.1094E-03 98.727 
31 2 .0000 .0000 2.4500 .02000 <-7.5316E-04 4.4391E-03 4,5025E-03 99.629 
32 2 .0000 .0000 2.4700 .02000 -8.8769E-04 4.7379E-03 4.8203E-03 100.612 
33 2 .0000 .0000 2.4900 .02000 -1.0256E-03 4.9554E-03 5.0604E-03 101.693 
34 2 .0000 .0000 2.5100 .02000 -1.1652E-03 “S.0891E-03 5.2Z06E-03 ~ 102.696 
35 2 .0000 .0000 2.5300 .02000 -1.3049E-03 5.1373E-03 5.3004E-03 104.252 
36 2 .0000 .0000 2.5500 .02000 -1.4430E-03 5.0989E-03 5.2992E-03 105.801 
37 2 .0000 .0000 2.5700 .02000 <-1.5778E-03 4.9737E-03 5.2180E-03 107.600 
38 Z .0000 .0000 2.5900 .02000 <-1.7076E-03 4.7623E-03 5.0592E-03 109.726 
39 2 .0000 .0000 2.6100 .02000 -1.8308E-03 4.4659E-03 4.8266E-03 112.292 
40 2 .0000 .0000 2.6300 .02000 -1.9460E-03 4.0864E-03 4,5261E-03 115.464 
41 2 .0000 .0000 2.6500 .02000 -2.0515E-03 3.6256E-03 4.1658E-03 119.502 
42 2 .0000 .0000 2.6700 .02000 -2.1461E-03 3.0856E-03 3.7585E-03 124.820 
43 2 .0000 .0000 2.6900 .02000 -2.2286E-03 2.4669E-03 3.3245E-03 132.094 
44 2 .0000 .0000 2.7100 .02000 -2.2979E-03 1.7679E-03 2.8993E-03 142.427 
45 2 .0000 .0000 2.7300 .02000 -2.3532E-03 9.8018E-04 2.5492E-03 157.387 
46 2 .0000 .0000 2.7500 .02000 -2.3937E-03 8.3523E-05 2Z.3952E-03 178.002 
47 2 .0000 .0000 2.7700 .02000 -2.4190E-03 -1.0547E-03 2.6389E-03 -156.441 
48 2 .0000 .0000 2.7900 -02000 -2.4287E-03 -2.7914E-03 3.7001E-03 -131.025 
49 2 .0000 .0000 2.8100 .02000 -2.4226E-03 -2.7923E-03 3.6968E-03 -130.945 
30 2 .0000 .0000 2.8300 -02000 -2.4010E-03 -1.0573E-03 2.6235E-03 -156.235 
31 2 .0000 .0000 2.8500 .02000 -2.3641E-03 7.9330E-05 2.3655E-03 178.078 
52 2 .0000 .0000 2.8700 .02000 -2.3124E-03 9.7436E-04 2.5093E-03 157.151 
33 2 .0000 .0000 2.8900 .02000 -2.2466E-03 1.7605E-03 2.8542E-03 141.917 
54 2 .0000 .0000 2.9100 .02000 -2.1675E-03 2.4580E-03 3.2772E-03 131.406 
35 2 .0000 .0000 2.9300 .02000 -2.0762E-03 3.0752E-03 3.7104E-03 124.025 
56 2 .0000 .0000 2.9500 .02000 <-1.9738E-03 3.6139E-03 4.1178E-03 118.642 
57 2 .0000 .0000 2.9700 .02000 -1.8617E-03 4.0734E-03 4.4786E-03 114.562 
38 2 .0000 .0000 2.9900 .02000 -1.7413E-03 4.4518E-03 4.7803E-03 111.362 
39 2 .0000 .0000 3.0100 .02000 -1.6141E-03 4.7472E-03 5.0141E-03 108.779 
60 2 .0000 .0000 3.0300 .02000 -1.4819E-03 4.9578E-03 5.1745E-03 106.641 
61 2 .0000 .0000 3.0500 .02000 -1.3461E-03 5.0823E-03 5.2575E-03 104.834 
62 2 .0000 .0000 3.0700 -02000 -1.2084E-03 5.1202E-03 5.2609E-03 103.280 
63 2 .0000 .0000 3.0900 .02000 -1.0707E-03 5.0718E-03 5.1836E-03 101.921 
64 2 .0000 .0000 3.1100 .02000 -9.3452E-04 4.9381E-03 5.0257E-03 100.716 
65 2 .0000 .0000 3.1300 .02000 -8.0147E-04 4.7209E-03 4.7884E-03 99.635 
66 2 .0000 .0000 3.1500 .02000 -6.7310E-04 4.4227E-03 4.4736E-03 98.654 
67 2 .0000 .0000 3.1700 .02000 -5.5087E-04 4.0463E-03 4.0836E-03 97.753 
68 2 .0000 .0000 3.1900 .02000 -4&.3609E-04 3.5946E-03 3.6210E-03 96.917 
69 2 .0000 .0000 3.2100 .02000 -3.2990E-04 3.0698E-03 3.0875E-03 96.134 
70 2 .0000 .0000 3.2300 .02000 -2.3317E-04 2.4712E-03 2.4822E-03 95.390 
71 2 .0000 -0000 3.2500 .02000 -1.4630E-04 1.7909E-03 1.7969E-03 94.670 
72 2 .0000 .0000 3.2700 .02000 -6.6339E-05 9.5517E-04 9.5747E-04 93.973 
- - = POWER BUDGET - - - 

INPUT POWER = 3.6409E-03 WATTS 

RADIATED POWER= 3.6409E-03 WATTS 

STRUCTURE LOSS= 0.0000E+00 WATTS 

NETWORK LOSS = 0.0000E+00 WATTS 

EFFICIENCY = 100.00 PERCENT 

**kkk*k DATA CARD NO. 7 EN OO 0 0 0 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 


0.00000E+00 0.00000E+00 


RUN TIME = 1. 952 
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